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PREFACE 



This small book written by the late Lev Artsimo- 
vich, a Full Member of the USSR Academy of Scien- 
ces, is the lecture course he delivered for the physi- 
cists interested in plasma physics. The book presents 
the fundamental information on the high-temper- 
ature plasma physics. These, at present largely well- 
established results, comprise, in effect, the sum of 
knowledge indispensable for any physicist with 
a wide enough sphere of interests. 

However, it should be clearly understood that 
the structure of the high-temperature physics has 
just been outlined and is still being developed. 
A major contribution to this field of modern physics 
has been made by the Department of Plasma Physics, 
headed by L. Artsimovich, which is a part of the 
I. Kurchatov Institute of Atomic Energy. L. Ar- 
tsimovich headed the research work with the toka- 
mak-type installations, adiabatic traps, rapid pin- 
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ches and other machines designed for harnessing 
the high-temperature fusion reaction. He was direc- 
tly involved in actual work in a number of these 
research projects and in development of the correct 
conceptual system and understanding of the plasma 
phenomena. This work, of course, could not fail to 
be reflected in the book; when discussing the expe- 
riments with tokamak-type installations and the 
adiabatic traps, L. Artsimovich described in more 
detail the most recent results available at that 
time. By now, some material in the book has 
become somewhat dated; for instance, plasmas 
generated in Tokamaks and adiabatic traps have 
now more impressive parameters and we have now 
a markedly better understanding of the processes 
in them. However, since the general conceptual 
system has not undergone any major changes we 
have deemed it unwise to alter the original author's 
presentation of ideas and, therefore, no significant 
changes have been made in the original text. 

Hence, the material presented in the book com- 
prises, in the opinion of L. Artsimovich, the fun- 
damental results obtained through many years of 
experimental and theoretical research in the high- 
temperature plasma physics. 

Prof. B. B. Kadomtsev 



1. 

DEFINITION OF PLASMA 

The term plasma is applied to the ionized gas all 
(or a considerable part) of the atoms in which have 
lost one or several of their electrons and converted 
into positive ions. This definition of plasma as a spe- 
cific state of matter is only a tentative one, it will 
be elaborated upon below. Generally, plasma can be 
assumed to be a mixture of three components — 
free electrons, positive ions and neutral atoms (or 
molecules). 

Plasma is the most widespread state of matter in 
the universe. The Sun and the stars can be regar- 
ded as enormous lumps of hot plasma. The Earth's 
atmosphere is surrounded by a plasma envelope 
known as the ionosphere. The circumterrestrial 
space outside the ionosphere contains peculiar plas- 
ma formations — the so-called radiation belts. In 
the terrestrial environment we encounter plasma in 
laboratories and technology in various gas dischar- 
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ges since any gas discharge (lightning, spark, arc, 
etc.) invariably generates plasma. 

The research in plasma physics was typically 
stimulated by potential practical applications of 
the results. At first, the scientists were interested in 
a plasma as a peculiar conductor of electric current 
and as a light source. Nowadays, we should consider 
the physical properties of plasma from another view- 
point and new aspects of plasma will be revealed 
to us. First, the plasma is the natural state of mat- 
ter heated up to a very high temperature; and, 
secondly, the plasma is a dynamic system subjected 
to the action of electromagnetic forces. The new 
approaches in the studies of the plasma behaviour 
stem directly from the large-scale technological 
problems involving physics. The major problems 
are the controlled nuclear fusion and the magneto- 
hydrodynamic conversion of heat energy into 
electric energy. In the near future plasma physics 
may play a significant role in accelerator techno- 
logy. 

The studies of plasma phenomena are of not only 
practical interest. The plasma is a material medium 
consisting of a group of particles which interact be- 
tween themselves via the simplest mechanism, name- 
ly, the Coulomb electrostatic forces. A physicist 
has to understand the mechanisms of various pro- 
cesses which occur in a plasma proceeding from 
the known microstructure of plasma. The basic 
theoretical concepts here are extremely clear. Here 
we remain invariably within the domain of the 
classical physics since the quantum effects do not 
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play any noticeable part in the typical plasma. 
Nevertheless, the theoretical analysis of the plasma 
effects is far from being completed and many gaps 
still remain in our understanding of the plasma. 

The major experimental efforts at present are 
aimed at developing the techniques for generating 
plasmas with high parameters, namely, high 
temperature and high density. Our subject of stu- 
dy — the high-temperature plasma — is produced 
by ourselves and we try to provide optimum condi- 
tions for the plasma existence when the plasma is 
in the quasineutral stable state. 

Before continuing further, let us elaborate upon 
the definition of a plasma. The electric forces which 
bind the opposite charges in a plasma provide for 
its quasineutrality, that is, an approximate equali- 
ty of -the concentration of electrons and ions. Any 
separation of charges when electrons shift with 
respect to ions must give rise to electric fields which 
tend to compensate the perturbation. These fields 
are the higher, the higher is the particle concentra- 
tion and for a high-density plasma they are very 
high, indeed. 

To estimate a field produced by a perturbation of 
plasma neutrality, let us assume that charges have 
been completely separated in a certain volume so 
that it contains only one type of charges. The elect- 
ric field in this volume satisfies the Poisson equation 

div E = 4np 

where p is the charge density. 

If the linear dimensions of the volume are of the 
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order of x and the concentration of the charged par- 
ticles in the plasma is n, we obtain 

div E ~ F/x ~ Anne 

and, hence, 

E « Annex 

The variation of the plasma potential in the volume 
of charge separation is V ?a Ex » Annex 2 . 

Let us consider the following example. Assume 
that the fully ionized plasma is produced from 
hydrogen which is, initially, at the normal tem- 
perature and under thepressureof linmHg. A cubic 
centimetre of such a plasma contains 7-10 16 ions 
and electrons each. In these circumstances E is about 
10 u V/cm. Therefore, if quasineutrality is sharply 
disturbed in a volume with a diameter of about 1 mm, 
the resulting electric field is more than 10 lu V/cm 
and the potential difference across this volume is 
of the order of 10 9 V. Such a separation of charges 
is, clearly, unfeasible. Even in a plasma of a much 
lower density a sharp disturbance of quasineutrality 
in such volumes will be immediately suppressed by 
the developing electric fields. The electric field 
will push out the charges of one sign from the volume 
with separated charges and draw the opposite char- 
ges into it. However, in a small enough volume o 
plasma, quasineutrality may be violated if the 
electric field generated by an excess of particles of 
one sign is too weak to have a significant effect on 
the motion of particles. For a given concentration 
and temperature of the plasma, it is described by 
a characteristic linear parameter 6: when x <^ 6, 
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then separation of charges within a volume with the 
linear size x does not produce a significant effect on 
the motion of particles; when x ^> 8, then concen- 
trations of the opposite charges in this volume are 
almost equal. 

The characteristic length 8 can be estimated in the 
following way. When the charges in the volume of 
the linear size 6 are completely separated, the 
potential energy of a charged particle is of the order 
of the energy kT of the thermal (heat) motion of 
particles; here T is the plasma temperature in abso- 
lute degrees *. Thus, we obtain the relationship 

U --= eV ~ AnneW ~ kT 

Hence we find 

6 « [kT/(/mne 2 )] v2 (1) 

The parameter 6 is derived also from the analysis of 
electric field screening in a plasma. Assume that 
a "test" point charge q has been introducad into 
a plasma. The potential generated by this charge 
at a small enough distance r from it is qlr. However, 
polarization of plasma due to the electric field of 
the charge alters the form of the potential function 
at large distances from the charge. 

* The thermal energy of plasma particles can 
also be expressed in electron-volts using 
the relationship kT = e6/300 eV = 1.6-10- 12 e eV 
which yields T = 11 6009 eV. The parameter 6 
is referred to as the plasma temperature 
expressed in electron-volts. The temperature 
denned in this way is identical to the parameter 
characterizing the thermal energy. 
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In a state of statistical equilibrium, the spatial 
distribution of electrons and ions in the vicinity of 
the test charge is given by the Boltzraann formula 
(n ~ exp (—U/kT), where U is the potential ener- 
gy. U has opposite signs for ions and electrons). 
This can readily be seen to result in screening of 
the electric field in this region. The concentration 
of particles with the charges opposite to q is higher 
in the vicinity of the test charge q, that is, for com- 
paratively high values of U/kT; this should result 
in a sharp drop in the electric field. Calculations 
using the Poisson equation and the Boltzmann dis- 
tribution law show that at large distances from the 
charge q the potential decreases exponentially and 
the electric field is high only in the sphere with the 
radius of the order of 6 around the charge q. Debye 
was the first to introduce the characteristic length 6 
in his study of strong electrolytes. Later, this concept 
was applied to plasma physics. The parameter 6 is 
referred to as the Debye radius or the Debye length. 
Inserting the numerical values of the constants into 
Eq. (1) for 6, we obtain 

6»7(77») 1/2 (2) 

Here T is the plasma temperature which is assumed, 
so far, to be the same for the electron and ion com- 
ponents of the plasma. While the Debye radius 
describes the spatial scale of the decompensation re- 
gions, the lifetime of these regions is found by di- 
viding 8 by the velocity of faster particles (electrons): 

t = 8/U e = (kT/4nne 2 ) U2 (m e /kT) i/2 = (m e /Anne 2 ) U2 
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The parameter 1/t has the dimensions of frequency 
and is identical to the frequency of the electrosta- 
tic plasma oscillations which occur in plasma when 
the electron groups are displaced from the equilibri- 
um distribution. This frequency 

<o = (4nne 2 /m e ) i/2i (3) 

is known as the plasma or Langmuir frequency. 

The higher the density of plasma the smaller the 
sizes of decompensation regions and their lifetimes. 
In a high-density cold plasma quasineutrality can 
be violated only within a fairly small region. In 
a low-density hot plasma the Debye length can be 
considerably larger than the size of the plasma. If 
this is the case, the motions of ions and electrons 
are independent and no automatic equalization of 
concentrations of the opposite charges occurs. 

We can now redefine the term plasma in the fol- 
lowing way making use of the concept of the Debye 
radius — the plasma is an aggregation of freely 
moving particles with opposite charges, that is, an 
ionized gas is called the plasma if the Debye length 
is small compared to the size of the volume occupied 
by the gas. This definition was put forward by 
Langmuir who was the pioneer of the plasma science. 

The two plasma parameters we have introduced — 
concentration and temperature — should be com- 
mented upon. 

1. Generally, the electron and ion concentrations 
need not be identical since a plasma can contain not 
only singly charged but also multiply charged ions. 



15 



If the concentration of ions with a unit charge is 
denoted by n 1 , the concentration of ions with the 
charge of 2 by n 2 , and so on, the electron concentra- 
tion n e is given by n x -f- 2n 2 -f .... However, below 
we shall mainly discuss the plasmas in which the 
election and ion concentrations are identical. In 
particular, this is the case for the hydrogen plasma. 
Normally, it is not too complicated to take into 
account the effects of the multiply charged ions on 
the basic processes in the plasma. 

2. The concept of the plasma temperature T is 
valid only if electrons and ions have the same kine- 
tic energies. Generally, the plasma must be charac- 
terized by at least two temperatures — the electron 
temperature T e and the ion temperature T t . The 
plasmas which are generated in a laboratory environ- 
ment or in technology have typically the tempera- 
ture T e which is usually much higher than T t . The 
difference between T e and T t is due to an enormous 
difference between the masses of electrons and ions. 
The external electric power sources that generate 
a plasma (in various gas discharges) transfer the 
energy to electrons since it is the electrons that are 
the current carriers. Owing to collisions with rapidly 
moving electrons, the ions increase their thermal 
energy. In such collisions the fraction of the kinetic 
energy of electrons which can be transferred to the 
ion is no more than imjmi, where m e and wi 4 are 
the masses of an electron and an ion. The mean 
fraction of the electron energy transferred to the ion 
in a collision is even lower. Since m e is much less 
than mi, the electron has to undergo very many 
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(thousands) collisions to transfer^completely to the 
ions the excess energy it has. Since the thermal ener- 
gy exchange between electrons and ions proceeds 
concurrently with the energy transfer from the elec- 
tric power sources to electrons and at the same time 
the energy leaks from the plasma via various heat- 
transfer mechanisms, typical result is a sharp tem- 
perature difference between electrons and ions in 
an electric discharge in gases. This energy difference 
decreases, typically, with increasing plasma concen- 
tration since the frequency of collisions between 
electrons and ions in a given plasma volume increa- 
ses with the squared concentration. 

Under some special conditions T t can become 
much higher than T e in a strongly ionized plasma. 
For instance, such conditions exist in short-duration 
high-power pulse discharges which result in devel- 
opment and cumulation of shock waves in the 
plasma. 

2. 

MOTION OF ELECTRONS AND IONS 

IN PLASMA IN THE ABSENCE 

OF EXTERNAL FIELDS 

The character of this motion is determined by the 
interaction between the particles. In the highly 
ionized plasma the principal type of interaction 
between the particles is the classical Rutherford 
scattering in the Coulomb field. We shall deal 
with three basic types of elementary scattering pro- 
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cesses, namely, scattering of electrons by ions, 
electrons by electrons and ions by ions. Other 
elementary processes either occur with photon emis- 
sion and their probability is relatively low or they 
involve also the neutral particles and their signi- 
ficance diminishes with increasing degree of ioniza- 
tion. The processes of the first kind are illustrated by 
emission of bremsstrahlung electrons in the electron- 
ion collisions, the processes of the second kind are 
illustrated by the ionization phenomena and atom 
excitation by electron impact and by the charge 
exchange processes occurring between ions and atoms. 
If we deal with a non-hydrogen plasma we have, 
generally, to take into account the interaction be- 
tween electrons and ions in various energy states. 
The intensity of radiation emitted by the excited 
ions can prove to be very high and can make a signi- 
ficant contribution to the energy balance of the 
plasma processes. We shall limit ourselves to ana- 
lyzing particle interactions in a fully ionized 
plasma. 

Let us assume that a charged "test" particle is 
travelling in the plasma (such a particle may be any 
plasma electron or ion whose path we shall trace). 
When travelling, this particle will be scattered in 
the Coulomb field of plasma electrons and ions which 
it encounters along its path. If we consider a light 
particle travelling among heavy ones (an electron 
among ions), then the scattering centres are assumed 
to be stationary. In this case the probability of 
scattering into a given angle is determined by the 
classical Rutherford formula. 
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Every scattering of a test particle travelling in 
the vicinity of the scattering centre results in the 
particle's trajectory being rotated by an angle 0, 
that is, its velocity along its initial direction of 
motion decreases from v to ycos 8. In most cases scat- 
tering occurs at large distances and typically results 
in a very small change in the trajectory (this is typi- 
cal of the Rutherford scattering in the electric 
field of point charges!). Therefore, here we cannot 
visualize, as is customarily done in the kinetic theo- 
ry of gases, the trajectory of a particle as being 
a broken line composed of straight segments or 
"paths" which connect the points of "collisions". 
Instead, we have a gently twisting line whose direc- 
tion varies owing to the effect of the numerous but 
very weak impetuses due to "collisions" with other 
particles. In effect, these impetuses generate a con- 
tinuous influenceof the plasma "microneld" on a tra- 
velling particle. This microneld is a superposition 
of the electric fields of individual particles. 

It is natural here to introduce the concept of the 
free path % of a particle defining it as the distance at 
which the particle changes its initial direction of 
motion. This definition yields the following equa- 
tion: 

dv = — vdx/k (4) 

Here do is the mean variation of the velocity com- 
ponent along the initial direction of motion when 
travelling the distance dx. This definition can be 
used to express % in terms of the integral over the 
angular distribution of scattered particles. If the 
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■velocity vector rotates by the angle alter colli- 
sion, then the velocity component along the initial 
direction of motion decreases by v (1 —cos 6). When 
a particle travels the distance dx in the volume 
with n scattering centres per 1 cm 3 , the particle 
is scattered nf (0) dQ dx timesmto the angle close to 
within the solid angle dQ. " The function / (0) is 
known as the effective differential^ cross-section of 
scattering and is determined by the nature of inter- 
action forces. Hence, the mean value of dv is 
given by 

dv = - vndx \ f (9) (1 - cos 6) dQ (5) 

From Eqs. (4) and (5) we obtain 

1/JU-n j /(G)(l— cosl))dQ. (6) 

If the scattering centres are assumed to be stationary 
point charges, then the function / (0) is given by 
the Rutherford formula 

/ (9) = 1/4 (flWs/mxD*) 2 (1/sin* 0/2) (7) 

where g t is the charge of the test particle, m^ is its 
mass and q 2 is the charge of the scattering centre. 
Substitution of Eq. (7) into Eq. (6) and integra- 
tion over angles from a minimum angle mln to it 
yield 

X =, i/Ann {m^lq^f (1/ln (2/0 mln )) (8) 

The value of m)n may be estimated as follows. The 
electric field of the scattering centre will be assu- 
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med to be a Coulomb one only at distances smaller 
than the Debye radius 6. At larger distances the 
potential decreases exponentially. Therefore, we 
have to ignore the collisions in which the particle 
passes the scattering centre at distances larger than 8. 
Now/ note that the distances of approach of the 
order of 8 correspond to extremely small scattering 
angles in the plasma. Under such conditions the 
trajectory of the test particle is almost a straight 
line and the closest approach distance is practical- 
ly identical to the impact parameter. For small- 
angle scattering in the Coulomb field the impact 
parameter b is related to the angle by the well- 
known equation 

=- qtfJrmPb (9) 

Substitution of b -- 8 into Eq. (9) yields 8 min . 
Since this calculation yields a very high log quantity 
in Eq. (8) for interactions of particles in the plasma 
(in alTthe cases of interest it varies from 10 4 to 10 8 ), 
a rough approximation of the estimated minimum 
scattering angle practically does not affect the com- 
putation of K. 

The assumptions used for deriving the free path X 
are valid when the test particle is an electron inter- 
acting with plasma ions. Denote by % Pi the mean 
free path of electrons colliding with plasma ions. 
The path % ei is derived by averaging Eq. (8) over 
the energy spectrum of electrons. When all the 
plasma ionsare singly charged an 1 the energy distri- 
bution of electrons is ffiven by the Maxwell formula, 
we obtain the following expression for the mean 



21 



free path: 

X e «^4.5.10 5 -f-f- (10) 
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Here L c is the so-called Coulomb logarithm which 
is derived by substituting into In 2/9 mln the mini- 
mum angle found from Eq.' (9) for b = 8, mv 2 = 
= 3kT e and q % = — q x = e. The value of L c 
varies from 10 to 20 for very wide variation ranges 
of n and T e . Since fairly rough estimates of the 
parameters of particle collisions can be used in 
plasma physics we shall assume below that L c =- 
== 15. Several otber parameters describing the col- 
lisions between electrons and ions can also be intro- 
duced in addition to X ei . The effective cross-section 
a ei for such collisions is given by the relationship 
X ei = l/rw ei , and the mean time period between 
two successive collisions is x ei — X e i/v e , where v e 
is the mean thermal velocity of electrons. The colli- 
sion frequency v el is the reciprocal of x ai . The above 
parameters may be calculated from 1he following 
expressions (for L c = 15): 

(11) 

The above parameters can be readily modified for 
describing the collisions with multiply charged ions. 
For these processes the effective cross-section o ei 
increases as the squared ion charge and the other 
parameters vary in a similar way. 

The collisions between electrons and ions arc the 
most, significant particle interactions in the plasma; 
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in particular, they determine the mechanisms of 
such processes as electric conduction an^ diffusion. 

To make a more comprehensive description of the 
Coulomb interaction between the plasma particles, 
we have to introduce also the parameters which 
will characterize the statistical effect of the colli- 
sions between identical particles (electron-elec- 
tron and ion-ion collisions). Here the calculations 
are more complicated since the analysis of the 
elementary collisions should take into account the 
motion of scattering centres. However, this can 
affect only the numerical factor in the mean free 
path expression, while the temperature dependence 
should remain the same. For instance, the expres- 
sion for X PK (the mean free path for the electron- 
electron collisions) must be the same as the expres- 
sion for % ei differing only by a numerical factor 
fairly close to unity. The expression for % tl (the 
mean free path for the ion-ion collisions) is derived 
from the formula for % ee by substituting T { for T e . 
The parameters x ee and x el are close to each other. 
The ratio x n h ee is Ym t /m e Vff/ff. When the 
electron and ion temperatures are equal, the ion- 
ion collisions occur much rarer than the electron- 
electron or electron-ion collisions. 

Let us summarize. The above analysis" has attem- 
pted to describe the interactions of charged parti- 
cles in the plasma in terms of the elementary kinetic 
theory of gases; we have replaced the smooth gradu- 
ally weaving trajectories of electrons and ions by 
arbitrary broken lines and reduced the statistical 
effect of numerous weak collisions to one abstract 
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strong impact. The usefulness of such, not too rigor- 
ous approach, lies in the fact that the expressions 
for the mean free path, the mean time between two 
successive collisions, etc. help to visualize the basic 
physical processes in the plasma. There exists 
quite a rigorous technique for analyzing the Coulomb 
interactions between the plasma particles using 
the mathematics of the kinetic equations theory. 
Let us now discuss the exchange of thermal energy 
between electrons and ions in the plasma. Let us, 
at first, consider the simplest case. Assume that 
a fast electron with the momentum p — m e v e 
travels in the vicinity of a stationary ion and is 
scattered into the angle 6. In the process the ion 
receives the momentum Ap =2p sin 0/2. As a re- 
sult, the ion starts moving with the kinetic energy 

AW = (2p;sin 9/2) 2 /2m; (12) 

To find the kinetic energy transferred by a fast 
electron to stationary ions per unit time, multiply 

AW by nvj (0)iQ and integrate over angles. The 
result is 

dW e fdt = Le'(4nne l fmtv,) (13) 

(it is assumed here that the ions are singly charged). 
The expression for the transferred energy can be 
rewritten as] 

Lc(4nne*/niiV e ) = (2m c fmi)'(4nne l /rrif,vf) (m,,y*/2) — 
= (2m ( ,/m i yv ei W e (14) 
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Here v ei is the number of collisions between the 
electrons with the kinetic energy W e and stationary 
ions per unit time. The mean relative fraction of 
energy lost in one collision is 2mjm i as should be 
expected from the model of the elastic collision 
between two spheres. The mean energy transferred 
from the plasma electron to the ions in 1 s is found 
from Eq. (13) by integrating over the Maxwell 
velocity distribution *: 

- dWjdt = 1.2- \Qr"n/AVTl (15) 

where A is the atomic mass of the ionized gas. This 
equation is valid only for T e ^> 7^; if T e is com- 
parable to T t, Eq. (15) must be replaced by the 
following equation: 

- dWjdt - (1.2 ■ lO-"n e /A) (T e - T t )ITl 12 (16) 

To estimate the efficiency of heat exchange between 
ions and electrons in the plasma, let us discuss 
a specific case. Let the electron temperature in the 
hydrogen plasma be maintained at 10 6 K and the 
electron concentration be n = 10 13 . Under these 
conditions heating of the plasma from K to T t — 
— ill T e takes 1.2 ms. Note that these calculations 
are meaningful only if in the time period under 
consideration the ions retain all the energy they have 
obtained. Since an ion with an energy corresponding 
to a few hundreds of thousand degrees travels about 
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100 m in 1 ms, the heating of ions under such condi- 
tions necessitates efficient thermal insulation of the 
plasma. 

A major characteristic of matter is the equation 
of state, that is, the relationship between the pres- 
sure, density and temperature. The equation of 
state for the plasma with isotropic distribution of 
velocities of charged particles has the same form as 
the equation of state for a two-component ideal 
gas: 

p = n e k (T e + T t ) (17) 

where p is the plasma pressure which is the sum of 
electron and ion pressures. The equation of state is 
known to have such a simple form only if the mutual 
potential energy of particles is negligibly small com- 
pared to their thermal energy. In the plasma, this is 
practically always the case. Indeed, the potential 
energy of a charged particle in the fluctuating 
plasma microfield should be of the order of e 2 /a, 
where a is the mean distance from a given particle 
to the next one (a ~ lfy^2n e ). It. can be readily 
seen that the ratio 

(e*/a)/(kT)~yn e ey(kT) 

is much smaller than unity except when we consider 
a plasma with the particle concentration of 10 18 and 
higher and the temperature not lower than 10 4 'K. 
The energy distribution of the gas particles'" is 
usually assumed to have the Maxwellian'~form. 
This assumption cannot be "automatically applied 
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to the plasma. The Maxwellian energy spectrum is 
established owing to collisions between the gas 
particles. 

The Maxwell distribution in a wide energy range 
can be established for a given group of particles 
with an arbitrary initial energy distribution, that is, 
the Maxwellian "tail" for the particles with energies 
much higher than kT can grow only during a certain 
time period in which the particles must undergo, 
in the average, several collisions with each other 
and these collisions must be between the identical 
particles. For instance, the energy distribution of 
electrons will become practically Maxwellian after 
a time period which is about 10 times the mean 
time r ee between two successive electron-electron 
collisions. For ions the Maxwellian distribution is 
established after about ten ion-ion collisions 
Thus, when T t is about T e "maxwellization" for 
ions proceeds much slower than for electrons. How- 
ever, for each plasma component, maxwellization 
due to collisions between the identical particles 
proceeds faster than the development of thermal 
equilibrium between both components. Therefore, 
both ions and electrons in the plasma have the 
Maxwell distributions but with two different tem- 
peratures. At the same time, if the charged plasma 
particles have a short lifetime, the energy distri- 
bution of ions in the range of energies much higher 
than kT can differ sharply from the Maxwellian 
distribution (the number of such fast ions will be 
much less than it should be according to the Maxwell 
formula). 
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3. 

PLASMA BEHAVIOUR 
IN ELECTRIC FIELDS 

An electric field or a pressure gradient in a plasma 
gives rise to directed particle fluxes (flows) — an 
electric current in the first case and a diffusion flux 
in the second case. When electric current is flowing 
in plasma, the ions can be assumed stationary. 
The current is carried by electrons. In the simplest 
case of constant current there should be established 
an equilibrium between the force of the electric 
field acting on electrons and the decelerating force 
due to collisions between electrons and ions. The 
latter force is given by the mean directed momentum 
lost by electrons in collisions with ions per unit 
time. An electron has v ei collisions in 1 s and in 
each collision it loses the momentum m e u, where u 
is the directed electron velocity. Hence, the dece- 
lerating force is m e uv ei and the equilibrium is es- 
tablished when 

eE = muvgi (18) 

The current density in plasma is 

/ = neu (19) 

Hence we obtain 

/ = ne*E/m c v ei = neH e iE/m e (20) 

Eq. (20) is the Ohm law for a plasma. The parameter 

Gx^nehei/me (21) 
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is the plasma conductivity. Insertion of the expres- 
sion for r ei and the numerical values of constants 
yields for L c = 15 

o E « 10 7 7f 2 CGSE units (22) 

This formula is valid for a fully ionized plasma 
with singly charged ions (the hydrogen plasma). 
Note that conductivity grows rapidly with increas- 
ing T e . At temperatures about 10 8 K, the conduc- 
tivity of the hydrogen plasma must be higher than 
the conductivity of copper at room temperature 
more than by an order of magnitude. If plasma 
contains multiply charged ions its conductivity is 
considerably lower. In the general case when the 
plasma contains ions with the charges z 1? z 2 , z 3 , . . . 
. . ., z n with the concentrations a l7 a 2 , a 3 , . . ., a„, 
Eq. (22) for the conductivity is transformed into 



2 a k z h 

0^1- 10 7 ^ -TT (23) 

h a * z l 
k 

Let us now discuss briefly the application of the 
Ohm's law to a plasma/Thislaw is valid if in plas- 
ma there is an equilibrium between the forces of the 
electric field acting on electrons and the decelerating 
forces. However, can such an equilibrium be estab- 
lished under any conditions? 

The decelerating force action on an electron acce- 
lerated by the field is the lower the higher is the 
electron's velocity. Take an electron from the far 
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tail of the Maxwell distribution (W e > kT e ). The 
directed velocity component acquired by the elect- 
ron between two "collisions" with ions is proportio- 
nal to x ci and, therefore, increases as y 3 . Therefore, 
if the thermal velocity v of our electron is high 
enough, its directed velocity u can become as high 
as v or even higher. Under these conditions it is no 
longer possible to use a simplified model in which 
the electron accelerates to a relativity low directed 
velocity along its free path and then completely 
loses it in an instantaneous strong "impact". Actu- 
ally, acceleration and deceleration of the electron 
are concurrent. While the electron is acquiring the 
directed velocity, the Rutherford scattering on ions 
is gradually changing the direction of its motion. 
The electric field tends to straighten the electron's 
trajectory while the interaction with ions is bend- 
ing it. If the increase in the directed velocity 
component is not compensated by scattering, the 
equilibrium of forces cannot be established and the 
electron should accelerate continuously so that its 
energy continuously increases. With increasing 
electron energy, the decelerating force decreases, 
hence, an electron involved in continuous accelera- 
tion by the field will be accelerated as long as it 
remains in the field. The field accelerates continu- 
ously those plasma electrons which have acquired 
along the free path X an additional velocity u which 
is higher than their initial velocity v. This condi- 
tion for continuous acceleration may be written as 
follows: 

eEx ei /m e >v (24) 



30 



Since x ei is proportional to tfln, Eq. (24) shows 
that continuous acceleration occurs only when 
EWJn is higher than a certain boundary value. As 
can easily be seen, for the hydrogen plasma this 
boundary value is of the order of 3-10" 12 , if E is 
measured in V/cm and W e in eV. In the plasma ex- 
periments the condition (24) is typically satisfied 
only for electrons with the energy much higher than 
kT e . The share of such electrons in the electron com- 
ponent is very low. In this case, the current com- 
prising the overwhelming majority of plasma elect- 
rons satisfies the Ohm law. However, there will 
exist in the plasma the current of the accelerated 
electrons for which the Ohm law is not valid. For 
large values of Eln, the condition (24) will be 
satisfied also for the electrons with the medium 
thermal velocity. In this case, continuous accelera- 
tion can involve the major part of the plasma elec- 
trons and the Ohm law will be sharply violated. 
As calculations show, the continuous acceleration 
process develops markedly when the ratio of the 
mean value of u for the plasma as a whole to the 
mean thermal velocity of electrons is greater than 0.1. 
The ratio ulv increases with v 2 so that when ulv 
for the electrons with the mean thermal velocity 
is 0.1, for the electrons with the energy of 10 kT e 
the velocity u is close to v. Such electrons are close 
to the threshold of continuous acceleration. It may be 
assumed that continuous acceleration of electrons 
is observed in ring-type electric discharges when 
a plasma is generated inside a toroidal chamber and 
accelerated by a cyclic electric field. In these 
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experiments under special conditions a certain 
(relatively small) group of plasma electrons has 
been observed to be accelerated to very high ener- 
gies at comparatively low voltages across the plas- 
ma ring. 

A more detailed analysis of the accelerated elect- 
ron flows indicates that these flows can generate 
and amplify a variety of oscillations and waves 
in the plasma and transfer their energy to them. 
This gives rise to a new mechanism of braking the 
accelerated praticles which puts a stop to accelera- 
tion after the plasma electrons have acquired a 
certain additional energy of a directed motion. 
This automatic braking mechanism makes it impos- 
sible for all the plasma electrons to be continuously 
accelerated. However, Eq. (22) cannot be used in 
this case for determining the plasma conductivity 
since braking of electrons through interaction with 
the waves must increase the resistance. The measu- 
red plasma conductivities in the ring systems agree 
with the above considerations. For high plasma den- 
sity and relatively low electric field the experimen- 
tal values of a E coincide to the measurement accu- 
racy with the values found from Eq. (22). The 
resistances are found to be abnormally high in low- 
density plasma at high temperatures. 

The peculiarities of the plasma are especially 
noticeable when we consider its behaviour under 
the effect of high-frequency electric fields. In this 
case the mechanical inertia gains particular signi- 
ficance. Let us consider the simplest case. Let the 
plasma be under the effect of the electric field 
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E = E exp (iat) (complex notation simplifies the 
analysis). If the held frequency is high enough that 
during one period of field variation the probability 
of collision between an electron and ions is small 
enough, we can neglect the braking force in the 
first approximation in our analysis of the motion 
of electrons. Hence the equation of motion can be 
written as 

m e x = — eE exp (ia>t) (25) 

where x is the coordinate along the field. Integration 
of Eq. (25) yields the velocity 

u ---= — x = eE exp (mt);'(im e a>) (20) 

Hence we obtain 

i = neu=- — E (27) 

Eq. (27) can be rewritten as 

E = im^j/ine 2 ) (28) 

Thus, the phase of a high-frequency field is shifted 
by 90° ahead with respect to the plasma current. This 
means that plasma in a high-frequency field posses- 
ses its own "non-magnetic" inductance which is due 
to the inertia of electrons. For low concentrations 
the non-magnetic inductance of a plasma conductor 
can be higher than its normal ("magnetic") inductan- 
ce. Let us analyze the conditions required for that. 
For a cylindrical uniform plasma conductor with 
the radius a, the non-magnetic inductance per unit 
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length is given by (in the GGSE units) 



E I nic m P 



J na 2 ne z Ne* 



where N is the number of electrons per 1 cm of plas- 
ma length. The normal inductance (self -inductance) 
for a conductor 1 cm long is of the order of unity 
in the GGSM units. In the GGSE units the induc- 
tance is of the order of 1/c 2 . Hence, the ratio be- 
tween the normal inductance and the non-magnetic 
inductance for plasma is 

iYe 2 /w e c 2 «3-10-'^ 

This ratio has a very simple physical meaning: it is 
the full number of electrons contained in a plasma 
conductor segment with the length of the classical 
electron radius (r e = e 2 /wigC 2 ). 

Eq. (27) relating the current density to the high- 
frequency field can be extended to take into account 
the braking of electrons due to collisions. Thus, the 
following relationship is obtained between E and /: 

E = / (p + mL e ) (29) 

Here p = II a E and L e = ro e /(rce 2 ) is the non-magne- 
tic inductance. 

Some comments should be made here on the die- 
lectric properties of plasma. These properties are 
ai o related to the motion of electrons under the 
effect of external electric field. In the above simple 
case of the external electric field E — E exp (icoi), 

the acceleration x and the velocity x of an electron 
are given by Eqs. (25) and (26). Further integration 
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over time yields the displacement 

x = eE exp (jco<)/(m e co 2 ) (30) 

Eq. (30) shows that the electron displacement x has 
a phase shift of 180° with respect to the force — eE 
acting upon the electron. This effect is the reverse 
of that encountered in normal solid dielectrics. 
When the phase shift between the displacement of 
the charge and the force acting upon it is 180°, the 
polarization of the material is opposite to the field 
and, hence, the dielectric constant is less than unity. 
The dielectric constant can be expressed in terms 
of the electric moment p of the unit volume in a well- 
known way: 

e = 1 + Anp/E (31) 

We have in the plasma 

p = — nex (32) 

where a: is the displacement of electrons (the displace- 
ment of ions can be neglected since it is very small). 
Eqs. (30)-(32) yield 

e = l-coo/(o 2 (33) 

Here co is the plasma frequency given by Eq. (3). 
When © is less than a> , the dielectric constant is 
negative. A consequence of this is that the electro- 
magnetic waves with the frequencies lower than <o 
cannot penetrate the plasma and are totally reflected 
from its surface. 

The diffusion processes in a plasma are, to a cer- 
tain extent, similar to those in normal gas. A diffu- 
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sion coefficient can be determined for each plasma 
component; similar to the kinetic gas theory the 
diffusion coefficient is ~kv t l'd, where A, is the mean 
free path and v t is the mean thermal velocity of 
the species. 

However, the plasma quasineutrality condition 
necessitates that the electron flux be equal to the 
ion flux across any surface in the plasma. Since 
electrons have a much higher diffusion coefficient 
than ions (the ratio is the same as the ratio be- 
tween the velocities for equal free paths), an electron 
concentration gradient at first will give rise to an 
electron flux from a higher-concentration region to 
a lower-concentration region greater than the ion 
flux. This will result in polarization of the plasma 
and in development of an electric field that will 
brake the electrons and accelerate the ions. Simple 
calculations show that the electron flux will be 
greatly decreased and the ion flux will be somewhat 
increased. This effect is known as the ambipolar 
diffusion. It is frequently found in the experiments 
with relatively weakly ionized cold plasma with 
diffusion of electrons and ions through neutral 
gas. 

The heat conductivity of plasma is determined by 
energy transfer in collisions between particles with 
different thermal energies in a region with tempera- 
ture gradient. The decisive part in heat transfer 
should be played by the collisions between the 
electrons since their collision frequency is high. 
The electron heat conductivity coefficient in a fully 
ionized plasma with singly charged ions can be cal- 
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culated from the following formula: 

S h «1.2-10-<T 5/2 (34) 

where 6/, is expressed in erg/cm • s • deg. Eq. (34) is 
valid when the electron mean free path is small as 
compared with the size of the region where the tem- 
perature gradient is established. Note that the rela- 
tionship between the heat conductivity coefficient 6^ 
and the electric conductivity coefficient a B is given 
by the well-known Wiedemann-Franz law. The heat 
conductivity of a fully ionized plasma increases very 
rapidly with temperature. Even when T is about 
10 5 K the heat conductivity of the hydrogen plasma 
is higher than the heat conductivity of silver at 
room temperature. 



PLASMA BEHAVIOUR 
IN MAGNETIC FIELD 

A physicist making his first acquaintance with the 
basic plasma concepts can readily feel some disap- 
pointment at the beginning since every thing seems to 
be too clear and almost self-evident. The work in 
plasma physics seems to be lacking in excitement, 
it looks like doing exercises with already known 
solutions. However, this attitude cannot develop 
when one considers the behaviour of a plasma in 
a magnetic field. 

Due to the effect of the magnetic field the plasma 
becomes non-isotropic and radically changes all its 
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properties. A magnetic field can be used to establish 
closed plasma configurations contained in a limited 
volume and, so to say, suspended in a vacuum. This 
is something that is quite untypical of gas. Under 
such conditions the plasma properties are closer to 
those of liquids. At the same time such a plasma 
possesses the properties which distinguish it from 
all other states of matter. Ultimately, all the 
characteristics of plasma processes depend on the mo- 
tion of particles; therefore, before we attempt the 
macroscopic description of magnetic properties of 
plasma, we have to consider in brief the effect of the 
magnetic field on the motion of electrons and 
ions. 

A charged particle in a uniform magnetic field is 
generally known to travel along a helical line. The 
projection of the trajectory on the plane perpendi- 
cular to the magnetic induction vector B is a circle 
with the radius p = mv L cl{qB), where v x is the 
transverse component of the particle's velocity. 
This is the so-called Larmor circle. The motion along 
this circle is the rotation with the Larmor frequen- 
cy co fl = qBI(mc). The particle" also travels along 
the lines of force with the constant velocity V||. 

Generally, the magnetic field is not uniform. 
However, in non-uniform fields typically considered 
in plasma physics, induction B has almost a con- 
stant value and sense at the distances of the order of 
Ihe Larmor radius of particles, that is, microscopic 
variations of the magnetic field are very small. Let 
us find how such weak non-uniformity of the field 
affects the motion of particles. 
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At first, assume that the field varies along the 
field line. The trajectory of the particle travelling 
along this line noticeably changes its shape within 
a segment in which the magnetic induction B in- 
creases or decreases significantly. When the particle 
travels towards the increasing field, the trajectory 
becomes more steep and it can be compared to 
a spring being compressed. When a particle travels 
towards a decreasing field, its trajectory becomes 
less steep. 

This effect can be readily explained. A charged 
particle which rotates along the Larmor circle 
gives rise to a circular current so that it is equiva- 
lent to an elementary diamagnetic with the magnetic 
moment ]x = W x /B, where W \ is the kinetic energy 
of the transverse motion. Indeed, the Ampere theory 
yields 

,. - l - n 2. I® [A- 1 QBm*v\c* _mv\_W ± 

M " _ c l " 2c P ~ 2c mcq*B*- ~ 2B ~ B { ° > 

The magnetic field whose intensity "varies along the 
lines of force acts upon a diamagnetic with the force 

F = - yidB/dl (36) 

(here differentiation is done in the direction of the 
field). The action of this force results in variation 
of the longitudinal velocity v^ according to the 

relationship: 
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Multiplication of both sides of Eq. (37) by v\\ yields 

dt ~ B dl dt ~ B dt ( ' 

When a particle travels in a magnetic field, W ± + 
+ W\i is constant. Therefore, Eq. (38) can be trans- 
formed into 

^± = ^™_ (39) 

dt B dt v ' 

Hence, we find 

dWJW = dB/B; WJB = const (40) 

Thus, the ratio W± IB is constant when a charged 
particle travels in the magnetic field whose intensity 
variation along the field lines is not too sharp. 

The constant W ± /B is usually referred to as the 
adiabatic invariant of motion. This emphasizes the 
fact that the particle travels in a slowly changing 
magnetic field. The kinetic energy of the transverse 
motion is W ± = W sin 2 a, where W is the full 
energy of the particle and a is the angle between 
the velocity direction and the field line. Since W 
is constant and WJB is the adiabatic invariant, 
the ratio sin 2 a/fi is also an adiabatic invariant. 
This shows that the slope increases with increasing B 
arid in the region with a stronger magnetic field the 
helical trajectory becomes more steep as discussed 
above. 

Assume that at a certain trajectory point a = a 
and B = B . Using these initial conditions, we 
can find the value of a at any trajectory point from 
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the relationship sin 2 oc/fi = sin 2 a /5 , that is, 

sin a =^ sin a "j/5/Z? (41) 

When a particle travelling towards the increasing 
field comes to the point where B = B /sm 2 a , 
then the angle a becomes 90° and hence the longitu- 
dinal velocity V\\ vanishes. This means that at this 
point the direction of the longitudinal motion is 
reversed. The particle is reflected from the region 
of high field towards the region of a lower field. 

Thus, high-field regions under certain conditions 
can act as some magnetic mirrors for charged parti- 
cles. For instance, if the field increases in opposite 
directions from a certain intermediate region, a 
charged particle can be blocked between two magnetic 
mirrors; it will oscillate along the field lines within 
a restricted space region. The particles with a large 
enough angle a (sin a > ]/ B min /B max ) will be 
confined in this region. 

Let us now consider the motion of particles in 
a non-uniform field whose intensity varies perpendi- 
cular to the field lines. Let us at first consider the 
simplest case when the velocity is perpendicular to 
the magnetic field. Fig. 1 illustrates the particle's 
trajectory under such circumstances. The magnetic 
field is perpendicular to the plane of Fig. 1. 

The field increases along the axis x. Here the par- 
ticle's trajectory projection in the xy plane is not 
a circle, since the Larmor radius at the right is less 
than the radius at the left. Clearly, the trajectory 
is not closed after a complete cycle of rotation. Each 
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Fig. 1. 

Drift of a charged particle in a non-uniform 

transverse magnetic field. 



rotation cycle of the particle comprises a loop and 
with each loop the particle spirals a distance Ay 
along the axis y, that is, perpendicular to the mag- 
netic field gradient. The particles of opposite charges 
spiral in opposite directions along the axis y. 
After several rotation cycles, it can readily be seen 
that the particle's trajectory is a path of loops 
stretching parallel to the axis y. Such a motion is 
known as the magnetic drift. The drift velocity of 
the particle is small compared to its Larmor velo- 
city (it is assumed here that variation of the field 
intensity at the distances of the order of the Larmor 
radius is small). 
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Let us note the following feature of the drift 
motion. In the course of this motion the particle 
does not travel towards stronger or weaker fields. On 
the contrary, the particle travels along a narrow 
path which stretches within a constant-field re- 
gion. Hence, the parameter WJB is an adiabatic 
invariant in the case of the magnetic drift. 

The magnetic drift velocity v d is ky/T, where T is 
the Larmor rotation period. The ratio between v d 
and the rotation velocity, clearly, must be of the 
order of p/l, where p is the Larmor radius and I 
is the characteristic dimension of non-uniformity 
of the magnetic field: 

Ml ~ (MB) (dB/dx) (42) 

Hence we find 

p dB mv\c dB //0 , 

Calculations which we omit here confirm Eq. (42) 
and yield the following expression for v d : 

1 mv\_c dB ,/.. 

Eq. (44) gives the velocity of the motion in the 
plane perpendicular to the induction B. However, 
the drift motion in the non-uniform magnetic field 
can also be related to the longitudinal velocity v\\ 
of the particle. The mechanism of such a drift is 
illustrated by Fig. 2. The thick lines in Fig. 2 
represent the field lines of the non-uniform magnetic 
field which are generally curved. The Lorentz force 
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Fig. 2. 

Drift of a charged particle in a non-uniform 
magnetic field with the velocity parallel 
to the field. 



is zero at the point M 1 where the particle's velocity 
is parallel to B. However, the particle will continue 
inertial motion and slip off the field line; this 
will give rise to a small velocity component perpen- 
dicular to the field at the point M 2 . The transverse 
velocity will automatically give rise to the Lorentz 
force. This force gives rise to a drift velocity which 
is perpendicular to the plane of Fig. 2. Considera- 
tion of special cases (for instance, the motion of 
particles in a magnetic field generated by a recti- 
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linear current) readily shows that the particle's 
drift due to the longitudinal velocity v\\ has the 
same direction as^he drift due to the transverse velo- 
city v L . 

The following expression can be derived for the 
drift velocity in the general case, when V|| ^= 
and v ± t^ 0: 

^ = ^■^11 + ^/2) (45) 

Here, (o c is the Larmor frequency and R is the radi- 
us of curvature of the field line. The velocity vector 
has the sense of the vector product [B X R]. Eq. (45) 
is valid when the following two conditions are met: 

v<i <€ V v \ + y 'fl» that is, the drift velocity is much 
smaller than the velocity of motion of the particle; 
and the current density in the region where the 
particle travels is zero (or low enough). When these 
conditions are satisfied, the relationship between the 
gradient of the field intensity and the radius of 
curvature of the field line is simple enough to make 
possible derivation of Eq. (45). If the second 
condition is not satisfied, the expression for the 
drift velocity has a more complicated form. Note 
that both components of the particle's velocity re- 
main constant in the course of the drift motion. Since 
the particle always drifts perpendicular to the gra- 
dient of /?, the value of B along the drift path re- 
mains constant. Therefore the drift motion does not 
violate the invariance of the parameler^VF^/jS. 

From the above it follows that the motion of 
a charged particle in the non-uniform magnetic 
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field can generally be expressed as a superposition 
of three types of motion: 

(i) rotation along the Larmor circle with the velo- 
city v x ; 

(ii) motion of the centre of the Larmor circle 
along the field line with the velocity i^; 

(iii) drift motion of the centre of the Larmor circle 
in the direction perpendicular to B and grad B. 

The instantaneous positions of the centre of the 
Larmor circle comprise the axial line of the trajec- 
tory which, at the same time, can be considered to 
be an averaged path of the particle. The form of this 
line is a fundamental geometric characteristic of 
the motion of a particle in the magnetic field. 

The nature itself provides us with an excellent 
illustration of the motion of charged particles in 
non-uniform magnetic fields, namely, the so-called 
earth radiation belts. The radiation belts were 
found by the first earth satellites and space rockets; 
they consist of high-energy electrons and ions 
trapped in the earth's magnetosphere. The motion of 
electrons and ions of the radiation belts in the mag- 
netic field of the earth is rather complicated. If we 
ignore the Larmor rotation and consider only the 
averaged trajectories, the resulting paths are shown 
by the thick line in Fig. 3. Charged particles in the 
magnetic field of the earth oscillate along the field 
lines and are reflected from higher-field regions near 
the magnetic poles. A slower drift due to the cur- 
vature of the field lines is superimposed on this 
oscillatory motion. The drift motion is in the azi- 
muthal direction. The particles of opposite charges 
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Fig. 3. 

Motion of charged particles in the magnetic 
field of the earth. The thick line 
is the averaged trajectory. 



travel in opposite directions going around the globe 
from east to west and from west to east. 

To end this brief summary of motion of particles 
in the magnetic fields, we should consider the situ- 
ation when, along with the magnetic force, there 
is a force of some other origin acting on the travel- 
ling particle. Of the highest interest is the case 
when such a non-magnetic force is perpendicular to B. 
We shall omit here very simple calculations which 
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show that a force F perpendicular lo B gives rise to 
the drift with the velocity 

^-^[(FXB)] (40) 

When the force is due to an electric field E, we 
have F = qE and 

^=^[EXB] (47) 

Here the drift velocity does not depend on the charge 
and sign of the particle. The drift motion of this 
type also does not violate the adiabatic invariance 
of W JB (the particle travels perpendicular to the 
force F and therefore its mean kinetic energy during 
one Larmor cycle remains constant). 

5. 

THE EFFECT OF THE MAGNETIC 

FIELD ON CHARACTERISTICS 

OF THE PLASMA 

Let us now analyze the effect of the magnetic field 
on the properties of plasma. Let us assume that 
a moderate-density plasma is in a strong magnetic 
held. Between two Coulomb collisions, each charged 
particle of the plasma travels in a helical trajectory 
along the field. If the held is uniform, the axial line 
of the trajectory is a field line. The particles can 
travel perpendicular to the field only owing to the 
Coulomb collisions. Each collision makes a particle 
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travel a distance of the order of the Larmor radius. 
If the collisions are rare (low-density high-tempera- 
ture plasma) the particles are, so to say, tied to the 
field lines. Such a plasma is known as a "magnetized" 
plasma. The degree of "magnetization" is given by 
the ratio kip, where k is the mean free path and p 
is the mean Larmor radius. When ^/p^> 1 (rare 
collisions and high field), the particle can travel 
an appreciable distance perpendicular to the field 
only if it travels a considerable distance along the 
field. When 'kip <^ 1 (high-density plasma, weak 
field) the motion of particles is practically isotro- 
pic; this means that the magnetic field has a weak 
effect on the behaviour of the plasma. 

The degree of magnetization can also be expressed 
in the following form: 

- , X qB 

r v mc D 

where co B is the Larmor frequency and t is the mean 
time between two successive collisions. The electron 
and ion components of the plasma have different 
magnetizations. Under normal conditions the value 
of (0 B x for the plasma electrons is much higher 
than for the plasma ions (when T t ~ T e , t he rati o 
between these parameters is of the order of ~\f m^m^. 
Hence, magnetization of electrons is higher than 
magnetization of ions. It can happen, for instance, 
that electrons are magnetized and therefore travel 
freely along the field lines, while the magnetic field 
by itself does not have a noticeable effect on the 
motion of ions. Under such circumstances, 
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electrons are tied to the magnetic field, while ions 
are confined to the same space region by the electric 
field of the electron plasma component. 

But we are primarily interested in high-tempera- 
ture fully ionized plasma; in such a plasma the ty- 
pical situation is when both electron and ion com- 
ponents are magnetized. The main parameters of 
experimental hot plasmas vary in the following 
ranges: concentration — from 10 10 to 10 14 , electron 
and ion temperatures — from 10 6 to 10 7 degrees, and 
magnetic field intensity — tens of kilogausses. For 
these ranges of n, T and B the magnetization para- 
meter (o b t for the electrons varies from 10 5 to 10 u 
and for the hydrogen ions, from 10 3 to 10 9 . 

Since a strong magnetic field restricts the motion 
of plasma particles in the plane perpendicular to 
the vector B, such a field can be used as a screen 
preventing the contact between the plasma and the 
walls of the vessel containing the plasma (Fig. 4). 
A cylindrical column of fully ionized plasma is 
inside a chamber with a strong magnetic field. The 
space between the plasma surface and the chamber 
walls contains only vacuum and the magnetic 
field lines. The electrons and ions of the plasma 
cannot penetrate this space. Thus, a strong magnetic 
field provides for an effective thermal insulation of 
plasma. But strictly speaking, this is not an equili- 
brium state. Owing to Coulomb collisions, the 
plasma will, sooner or later, spread all over the 
chamber's volume up to its walls. The lifetime of 
an insulated plasma column is determined by the 
diffusion rate of the plasma particles in the magnetic 
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Fig. 4. 

A plasma column in a magnetic field. 



field (in the plane perpendicular to B). The diffusion 
theory shows that this lifetime is of the order of 
a 2 /cr£) x , where a is the radius of the plasma column 
and Z> ± is the transverse diffusion coefficient. A very 
rough estimate of D x can result from the following 
simple calculations. During the period t of one 
free path, the particle in the average makes one 
Coulomb collision which results in the particle's 
displacement in the plane perpendicular to B to 
a distance of the order of the Larmor radius p. To 
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find the resulting displacement after numerous 
collisions, we have to add up the squares of the 
individual displacements according to statistical 
relationships. Thus, the displacement of a particle 
perpendicular to B during the time t is 

Az « p YTR (48) 

At the same time, the mean displacement during 
this diffusion process is of the order of YD L t. 
I fence, we obtain 

D ± ~ p 2 /x ~ h>/(<oT) 2 (49) 

The mean time between two successive collisions is 
proportional to T 3 * 2 and therefore we find 

D x ttAn/B 2 YT (50) 

For high values of B and T, the coefficient D x must 
be very small. But the motion of particles along the 
1 1 eld lines is be same as when 5 = 0. Thus, the 
plasma in a high magnetic field has a sbirply 
marked anisotropy with respect to diffusion. 

The above qualitative discussion fails to account 
for one important feature of the diffusion mechanism. 
Rigorous theoretical analysis shows that the trans- 
verso diffusion of plasma particles is due only to 
collisions between different particle species, that 
is, to collisions between ions and electrons. The 
collisions between the same species cannot give rise 
to a macroscopic variation of the concentration 
profile. 
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Rigorous analysis yields the following expression 
for the diffusion coefficient: 

Dx^^y— [ —^ = D or — i— (51) 

Here v e is the mean thermal velocity of the plasma 
electrons, u> Be is the electron Larmor frequency, 
and Z) is the diffusion coefficient in the absence of 
the field. 

Clearly, the heat conductivity of plasma in the 
direction perpendicular to B also must sharply 
decrease with increasing field intensity. In contrast 
to diffusion which is determined by the collisions 
between ions and electrons, the heat transfer across 
the field lines in plasma occurs primarily in the 
ion-ion collisions (if T t is not too small compared 
to T e ). This is explained by the fact that the inten- 
sity of heat transfer is determined by the size of 
the region in which there meet the trajectories of 
particles with different thermal energies under the 
conditions of the temperature gradient. 

The thermal conductivity coefficient in the direc- 
tion perpendicular to B is proportional to the squa- 
red size of this region and this size is of the order of 
the Larmor radius. Therefore, the ion component is 
primarily responsible for heat transfer. 

The coefficient of ion thermal conductivity in the 
direction perpendicular to the field lines decreases 
by a factor of the order of (toBita) 2 compared to 
the coefficient in the absence of the magnetic field. 
The coefficient k x of the transverse thermal conduc- 
tivity for the hydrogen plasma is given by the fol- 
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lowing expression: 

K 1 =2-10- 18 « 2 /(5 2 /T7) 

The rate of equalization of temperature along the 
radius in the cylindrical plasma column proves to 
be much higher than the rate of equalization of 
concentration (approximately, by a factor of 
Ymilm e ). Hence, a temperature gradient in the 
direction perpendicular to B in the plasma column 
must vanish long before the plasma spreads all 
over the space owing to transverse diffusion. 

Low diffusion rates make it possible, in principle, 
to establish in the magnetic field plasma configura- 
tions which are confined to a finite volume of space 
surrounded by vacuum on all sides and have long 
lifetimes (of the order of seconds or even tens of 
seconds). 

Let us attempt a macroscopic approach to analyze 
the behaviour of such confined plasma configurations. 
The first problem that arises here is as follows. The 
plasma has a gas-kinetic pressure p = nk (T e + 
+ Ti)\ therefore, a plasma configuration with 
a certain spatial distribution of pressure can exist 
for a long time only if in any plasma element the 
force due to the gradient of gas-kinetic pressure is 
balanced by forces of other nature. What are these 
forces? For instance, if we have a single plasma 
column with a constant pressure p, what is the 
force that compensates the pressure at the boundary 
of the column? 

In the light of the above discussion of the motion 
of charged particles in magnetic fields, the answers 
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to these questions seem almost self-evident. A con- 
fined plasma configuration must be contained by 
electrodynamic forces which arise due to the dia- 
magnetic behaviour of the plasma in the magnetic 
field. The plasma is diamagnetic since the Larmor 
currents of the rotating charged particles generate 
in each plasma element a magnetic moment directed 
opposite to the external field. This decreases the 
magnetic field in plasma and gives rise to a ponde- 
romotive force which compensates the pressure 
difference between the interior of the plasma and 
its boundary. 

The ponderomotive force action on unit volume 
of the plasma can be expressed in terms of the 
density / of the diamagnetic current in the plasma 
and the magnetic induction B. This force is [j X B]/c. 
The equilibrium condition has the following form: 

grad p = [jx Bl/c (52) 

Note that the current density j can be found as the 
sum of the elementary currents due to rotation of 
electrons and ions in the magnetic field B. 

The mechanism of plasma confinement can be 
illustrated by the following simple example. Fig. 5 
shows a diagram of the Larmor currents in a cylin- 
drical plasma column inside which Ihe pressure p 
is constant. Under these circumstances, the Larmor 
currents inside the plasma column cancel out each 
other so that / = 0. 

But at the boundary of the column there is a net 
circular current. The force resulting from interac- 
tion of this current with the magnetic field compen- 
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Fig. 5. 

Larmor currents in the plasma. 



sates the pressure difference at the boundary of the 
plasma column. 

Eq. (52) derived for equilibrium conditions can be 
extended to analyze the dynamics of plasma configu- 
rations. The second Newton's law for the motion 
of a unit plasma volume can be written as 

pdx/dt = [j x B]/c - grad p (53) 

where p is the plasma density and v is the velocity 
of the plasma volume. The right-hand member of the 
equation is the force acting on unit plasma volume. 
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Eq. (53) lacks all characteristics of the plasma as 
an ionized gas. This equation is generally known as 
the fundamental equation of magnetic hydrodyna- 
mics. It can be used for analyzing the processes 
which occur in any conducting medium in a strong 
magnetic field when the medium is capable of 
changing its form under the effect of external for- 
ces, for instance, for analyzing the phenomena in 
a high-conductivity fluid. In plasma physics 
Eqs. (52) and (53) give an approximate model of 
the physical process which is rather rough. How- 
ever, this model yields very'good results in the ana- 
lysis of statics and dynamics of plasma configura- 
tions if we ignore some finer features of the processes 
and some aspects of the dissipation phenomena in 
the plasma. 

The next approximation is given by a two-fluid 
model of plasma which accounts for the electron 
and ion components of the plasma. The following 
equations describe the dynamics of plasma proces- 
ses in the two-fluid model: 

Pe^- = -i*. (E + -^[v e x Bl) -giadp.-F,, (54) 
p i ^ = » e (E + 4-[viXB])-gradp 1 + F el (55) 

where p e and p ( are the densities of the electron and 
ion components, v e and Vj are the averaged veloci- 
ties of electrons and ions, E is the electric field in- 
tensity, and F CJ - is the friction force acting on the 
electrons in unit volume due to their collisions 
with ions. This force is given by the momentum 
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transferred per unit time from electrons to ions: 
F>t= P ,(v,.-v,) (5(}) 



1 n 



The physical meaning of Eqs. (54) and (55) is 
quite clear as the right-hand members of these 
equations are the net forces acting on the given 
particle species in unit volume of the plasma. 

The two-fluid model of plasma can be used to 
analyze such processes as electric current in a mag- 
netized plasma and the development of a variety 
of non-stationary plasma phenomena. It should be 
emphasized that both Eq. (53) and Eqs. (54) and (55) 
are valid only when the plasma pressure is isotro- 
pic so that p is a scalar quantity. This condition is 
satisfied only for such plasma formations in which 
plasma can freely spread along the lines of force of 
the magnetic field. 

6. 

RESULTS DERIVED FROM 

MAGNETOHYDRODYNAMIC 

EQUATIONS 

Let us now consider some results which can be 
derived from general magnetohydrodynamic equa- 
tions. A simple case to analyze is the equilibrium 
of a plasma column in the longitudinal magnetic 
field (see Fig. 4). The field lines here are straight 
and parallel to each other so that the field 
intensity has only one component, which we denote 
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as B z . The intensity B z is a function of the coordi- 
nates x and y which are perpendicular to the field. 
The equilibrium condition f or the plasma column is 

grad p = -1 [j x B] = ~ [rot B X B] (57) 

Hence we obtain 
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so that 

p + BVSn = const (59) 

Eq. (59) shows that the magnetic pressure B 2 /8n 
outside the plasma is higher than inside it by p. 
The maximum pressure at which the plasma can be 
confined by the field of the given intensity B is 
found from the condition 

P max =Bl/8n (60) 

In this case the plasma must completely push out 
the field from the volume occupied by the plasma. 
If the field lines are curved, the sum of gas-kinetic 
and magnetic pressures is not generally conserved. 
However, Eq. (60) regarded as the boundary con- 
dition remains valid. If the pressure at the curved 
surface of the plasma configuration drops to zero, 
then the following condition must be satisfied: 

p + Bf/8n = B 2 /8n 

Here B t and B are the field intensities at both sides 
of the boundary. 
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A high longitudinal current'passing through a plas- 
ma can produce a quasi-stationary cylindrical 
plasma column. In this case the plasma pressure is 
contained by the magnetic field whose circular lines 
surround the plasma column. Let us limit ourselves 
here to the simplest case when current flows only 
along the thin surface layer [of the plasma (strong 
skin effect). In this case the field intensity inside the 
plasma is zero and the pressure in it is constant. 
This pressure can be found from Eq. (60), where B 
is the intensity of the field generated by the current 
at the outer surface of the plasma. If the current is / 
(in the CGSE units) and the plasma radius is a, 
Eq. (60) yields 

1 / 2/ \ 2 _ /» 
P ^ &n [ ca / := 2jwV ^ ' 

Further, let us assume that the ion and electron 
temperatures in the plasma column are the same (T). 
We have here 

p = 2nkT (62) 

Eqs. (61) and (62) yield the following relationship 
between the temperature and the current: 

T = P/Ac 2 NkT (63) 

Here N — zinc? is the number of particles of the 
same species per 1 cm of the length of the plasma 
column. Eq. (63) has been derived with a special 
assumption that the current flows in a thin surface 
layer. But a more careful analysis shows that this 
formula is valid for any current distribution over 
the cross-section of the plasma column. 



GO 



The formation of a plasma column stabilized by 
the magnetic field of the current flowing in it is 
known as the linear pinch effect. A current passing 
through a straight discharge tube filled with gas 
should, at first, give rise to a plasma (owing to 
ionization of the gas) which will be compressed by 
electrodynamic forces so that a plasma column can 
be formed. In these circumstances the current does 
three things at once — it generates a plasma, heats 
it by the Joule heat and compensates its pressure in 
the compressed plasma column by the pressure of 
its own magnetic field. At the early stages of high- 
temperature plasma physics, this technique for 
generating high-temperature plasmas seemed to be 
highly promising owing to its apparent simplicity. 
It was thought that a sufficiently high pulsed current 
passing through a gas at low pressure can produce in 
a discharge tube over a negligibly short period of 
time a plasma column with a tremendous temperature. 

Indeed, assume that a current of 10 6 A passes in 
a discharge tube with a 10 cm diameter filled with 
hydrogen at the initial pressure of 0.01 mm Hg 
(the initial concentration of atoms n = 7-10 15 cm- 3 ). 
If all the gas is ionized and a plasma column is 
formed whose pressure is compensated by electro- 
dynamic forces, the plasma temperature, according 
to Eq. (63), must be 

T — A — ~ 3 . 1 fl'K 
400 Nk ~ 400 na*n k 

In practice, however, the highest temperature obta- 
ined with this technique was only about 10 6 degrees 
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(and that was only i'or a few microseconds). Experi- 
ments have shown that a high-current short-dura- 
tion pulsed discharge does not generate a quasi- 
stationary state described by the equilibrium condi- 
tion (57). Now it can be seen that we could have 
expected that if we had a keener insight. 

At the initial stage of discharge, the current flows 
only in a comparatively thin layer of the ionized 
gas near the walls of the discharge tube owing to 
the skin effect. Under these circumstances, the 
eiectrodynamic forces are not compensated by the 
gas-kinelic pressure of the plasma, since the cylin- 
drical layer with the current surrounds the practic- 
ally unperturbed weakly ionized gas. Then the 
plasma layer will be accelerated by eiectrodynamic 
forces towards the axis of the tube. This stage of 
the process can be described very roughly by 
Eq. (53) without the term grad p. The accelerated 
motion of the plasma towards the axis of the tube 
can be regarded as the development of a cylindrical 
compression shock wave. As the plasma layer travels 
towards the axis, its mass steadily increases and 
new gas layers before the plasma front become in- 
volved in the compression process. At the same time 
the gas becomes ionized. At the last stage of com- 
pression, the plasma accelerated by the magnetic 
pressure reaches the axis of the tube. At this stage, 
a considerable part of the kinetic energy of plasma 
converts into heat. Now it seems to be the time 
for the next stage of the process at which the bal- 
ance can finally be established between the electro- 
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Fig. 6. 

Deformation of a plasma column 
carrying current. 



dynamic forces and the internal pressure of the 
hot plasma. However, here comes into play anoth- 
er physical factor, namely, the instability of the 
plasma pinch. Under typical experimental condi- 
tions deformations of the "bottleneck" and "snake" 
types develop in the plasma pinch in a few micro- 
seconds (see Fig. 6). These deformations disrupt the 
regular geometrical structure of the plasma so that 
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Fig. 7. 

Compression of plasma by 

a rapidly increasing magnetic field. 



it starts to interact strongly with the walls of the 
discharge tube and rapidly cools down. 

This case illustrates for us for the first time a ma- 
jor problem of modern plasma physics, namely, the 
stability of plasma configurations. The plasma is 
very unstable in its behaviour; it tries in all possi- 
ble ways to get rid of the magnetic confinement 
slipping through the lines of the surrounding magne- 
tic field. At the end of the book we shall discuss the 
general factors which determine stability or insta- 
bility of plasma configurations. 

A rapidly increasing longitudinal magnetic field 
compressing the plasma provides for more effective 
heating of it than in the case of a linear pinch 
(see Fig. 7). This process also starts with the initial 
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non-stationary compression stage at which a shock 
wave travelling to the axis is generated. But when 
the plasma cylinder radius ceases to oscillate, 
a balance of forces is established for some time and 
adiabatic heating of plasma occurs with further 
increase of magnetic pressure. 

A magnetic field which increases up to about 
100 kOe in a few microseconds makes it possible to 
produce a short-lifetime plasma with a temperature 
about 10 7 degrees and a concentration over 10 16 cm -3 . 
The longer the discharge tube and the coil which 
generates the field the longer the lifetime of the 
high-temperature plasma. But as the plasma can 
flow out of the open ends along the field lines, the 
lifetime of high-temperature plasma is still short — 
for the above values of n and T it cannot be more 
than a few tens of microseconds in a system about 
1 m long. 



DEVELOPMENT OF EXPERIMENTAL 
TARGET 

The fundamental problem in the experimental 
physics of high-temperature plasma is to generate 
a medium to be studied. Serious experimental studies 
of the main plasma phenomena can be feasible 
only when there have been developed the techni- 
ques for generating high-temperature plasma whose 
quasistationary state can exist for relatively long 
periods. In recent years, some advances have been 
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made in development of such techniques. There 
are two main fields of research here. The first one 
is concerned with magnetic systems which make 
it possible to generate closed plasma configurations 
of a toroidal shape. Such machines are often called 
the toroidal magnetic traps. The second one is 
concerned with the systems in which the high-tem- 
perature plasma is confined in open mirror-type 
magnetic traps. In such traps the magnetic field 
increases along the field lines in both directions 
away from the plasma region. As discussed above 
the higher-field regions act as mirrors which reflect 
ions and electrons of the plasma. 

Let us discuss at first closed plasma configurations. 
In such a configuration the plasma spreads freely 
along the field lines. This means that the field lines 
must lie within the closed surfaces which are inside 
the chamber in the region where plasma is generated. 
A gas-kinetic pressure p of the plasma is a scalar. 
When the forces acting on the plasma are balanced, 
the pressure must be constant along the field line 
since as shown by Eq. (52) the pressure gradient is 
perpendicular to B. Superficially, there would seem 
to be a very simple technique for magnetic con- 
finement of the circular plasma pinch. This technique 
consists in winding a regular coil on the surface of 
a round toroidal chamber which generates a field 
with circular field lines (Fig. 8). This is the field 
that will confine the plasma which will be generated 
inside the chamber in a certain way. I 

However, this technique is essentially defective 
as can readily be shown by a simple theoretical argu- 
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Fig. 8 

Plasma ring in the field 

of a toroidal coil. 



merit. The magnetic field in this case is not uni- 
form, since B ~ HR so that ions and electrons of 
plasma will drift perpendicular to the field lines 
to the chamber walls. The balance equation (52) 
also shows that plasma cannot be confined in a sim- 
ple toroidal field. To support this assertion, let us 
discuss a special case. 

Let the plasma completely push out the magnetic 
field from the region it occupies. Then the following 
condition must be satisfied everywhere at the boun- 
dary of the plasma ring] — p — B 2 /8n. But the 
field is different at the points M x and M 2 (see Fig. 8) 
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Fig. 9. 

A figure 8 configuration. 



so that balance is impossible. It can be seen from the 
balance Eq. (52) that this is so in the general case 
when p =£/? 2 /8ji. The above simple technique 
for confinement of a circular plasma pinch can be 
modified in several ways. One modification is to 
make a figure 8 out of a toroidal chamber with the 
magnetic coil on it as shown in Fig. 9. Then the 
plasma pinch also has the figure-8 shape and the 
drift of particles has opposite directions at the 
opposite curved parts of the pinch, so that the 
displacement due to drift during one cycle will be 
cancelled out. Thus, the laws of motion of the 
plasma particles in the toroidal "figure 8" do not 
preclude equilibrium of the plasma. The results 
of the macroscopic analysis of equilibrium using 
Eq. (52) should be expected to show that the systems 
of this type can be utilized as magnetic traps for 
closed plasma rings. 

The toroidal figure 8 was, in fact, the first step 
in the development of a wide variety of magnetic 
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Fig. 10. 

Superposition of the fields Bq and B, v 



systems for generating closed plasma configurations. 
These systems based on an elegant concept ol 
Spitzer are known as stellarators. In stellarators 
the equilibrium of high-temperature plasma is 
maintained by using external magnetic fields. 

The stellarator studies have been for many years 
an important field of the high-temperature plasma 
physics. However, we cannot discuss here this very 
interesting work because this would lead us away 
from plasma physics towards the analysis of geo- 
metry of complex magnetic fields. 

Equilibrium toroidal configuration of another 
type can be generated using superposition of two 
fields (Fig. 10). One of the fields is produced by 
a toroidal coil (# ) and the other field B v is produ- 
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ced by longitudinal circular current flowing in the 
plasma itself. The field in this system has simple 
helical lines which wind about the axis of the 
plasma pinch where the magnetic field of the cur- 
rent vanishes. The drift displacements of particles 
are compensated owing to the curvature of the 
field lines. To understand the mechanism of this 
compensation without calculations, we shall illustrate 
it by a special case when the following two con- 
ditions are satisfied: (1) the magnetic field generated 
by the current in the plasma is low compared with 
the toroidal field generated by the external coils; 
and (2) a plasma particle has a sufficiently high 
longitudinal velocity and therefore can freely travel 
along the field line without being reflected from 
the higher-field regions. Under these circumstances 
the motion of the centre of the Larmor circle is the 
superposition of two motions — the motion along 
the field line (always in the same direction) and 
the drift due to non-uniformity of the toroidal 
magnetic field. Fig. 11 illustrates a simple model 
for superposition of these two motions in which 
the continuous drift along the axis y (the principal 
axis of the toroidal system) is expressed by two 
jump-like displacements. We see that the motion 
of the particle alorg the winding field line results 
in closing of the projection of its trajectory on the 
cross section of the plasma pinch so that the net 
drift displacement is zero. 

Analysis of the motion of particles in the super- 
imposed fields of the above type shows that drift 
displacements are compensated not only for the 
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Fig. 11. 

Model of the drift motion of particles 
in a helical field. 



particles with high longitudinal velocities which 
freely travel along the helical field lines, (the 
so-called "through-going" particles) but also for 
the so-called "confined" particles which oscillate 
along the field lines owing to the reflection from the 
higher-field regions. 

It should be stressed here that the compensation 
of drift displacement does not preclude any effect 
of the drift motion on the shape of the trajectory. 
Even the schematic model of Fig. 11 shows that 
the drift motion results in some deformation of 
the trajectory. This deformation is even more 
strongly marked for the trajectories of the confined 
particles. We shall discuss some of the resulting 
effects below when mentioning the problems of 
the rates of diffusion and heal transfer. 
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8. 

CONDITIONS FOR 
EXISTENCE OF A PLASMA RING 

Let us now discuss in more detail the conditions 
for the existence of a plasma ring in a longitudinal 
external magnetic field and the magnetic field gene- 
rated by the current flowing in the ring. We illustrate 
the general problem by analyzing in detail 
this special case because, first, this is the simplest 
example of a magnetic confinement and, second, 
the longest times of high-temperature plasma con- 
finement have been obtained in the experiments 
with the machines utilizing this concept. 

Before we analyze the forces acting on the plasma 
ring, let us look at some of the characteristics of 
the magnetic field geometry in the toroidal trap 
under consideration. As mentioned above, the 
helical field lines wind about the circular axis of 
the plasma ring. Let the plane S be the cross sec- 
tion of the trap. The field line through the point M x 
(Fig. 12) in this plane after passing around the 
toroidal magnetic system will intersect the plane S 
at the point M 2 , next time it will intersect it at the 
point M 3 , etc. The set of points M l5 M 2 , M 3 , . . . 
is, generally, infinite. But for some field lines this 
set can be finite so that they are closed by them- 
selves. Such lines are known as "degenerate" lines. 
The axis of the ring at which the field is zero is 
such a degenerate line; it is shown by point in 
Fig. 12. 
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The behaviour of the field line which winds re- 
peatedly about the toroid is determined by the 
positions of the "descriptive" points M ly M 2 , M 3 , .... 
Assume that the point M n after n cycles is close 
to the initial point M v Then after the next cycle 
the angle of rotation about O will be more than 
360°. Where will be the next point? It is natural 
to suggest that the point M n+1 will be between the 
points Afj and M 2 , the point M n+i between the 
descriptive points M 2 and M s and so on. Thus, after 
numerous cycles, the points in the cross section 
of the trap will make up a smooth closed curve 
(curve F in Fig. 12). Thus, we can say here that there 
is a toroidal magnetic surface which, so to say, con- 
sists of (or is tightly woven of) one infinite field 
line. 

Analysis of the equations for a magnetic field 
supports the assumption about magnetic surfaces 
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Fig. 13. 

Magnetic surfaces in a helical field. 



in the system with the external toroidal field B$ and 
the field B v of the plasma current. As shown by 
the analytical results, for regular distribution of 
the current density over the cross section of the 
plasma ring the sufficient general condition for 
the existence of magnetic surfaces is the axial 
symmetry of the field. When the axial symmetry 
exists, the equations for the magnetic field lines 
can be integrated yielding the analytical expres- 
sion for magnetic surfaces. The magnetic surfaces 
comprise the set of embedded toroids (see Fig. 13). 
Any toroid out of this continuum is comprised of one 
field line. However, this set of magnetic surfaces 
includes a countable set of "degenerate" surfaces 
generated by closed field lines. In this case the 
surface is woven of a continuous set of field lines 
shifted with respect to each other. 
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Spitzer has suggested that magnetic surfaces 
exist also for such toroidal fields in which the 
helical field lines winding about, the circular axis 
are generated by the external current sources 
using special helical windings or by deforming the 
whole magnetic system into, for instance, a figure 8. 

However, theoretical analysis shows that for 
such asymmetrical fields the concept of embedded 
magnetic surfaces can be used only to some appro- 
ximation. The structure of the magnetic field is, 
generally, rather complicated. The individual mag- 
netic surfaces are separated by a system of toroidal 
tubes ("filaments") between which the behaviour of 
the field lines is quite chaotic. Clearly, the existence 
of magnetic surfaces is very important since it 
determines the prospects for utilization of magne- 
tic systems of various types for confinement of hot 
plasma. Since the plasma freely spreads along the 
field lines, its pressure p must be the same at diffe- 
rent points of the same magnetic surface. Thus, 
magnetic surfaces are families of plasma isobars. 
If the points 71/ lt M 2 , M 3 , ... fill up a certain 
wide region in the cross section of the plasma ring 
instead of comprising a smooth curve, then the 
plasma within this region can be at equilibrium 
only if grad p = 0. The fact that the pressure 
gradient vanishes at the cross section of the plasma 
ring is, clearly, equivalent to the development of 
abnormally high diffusion in the direction perpen- 
dicular to B. 

Let us now return to the analysis of the equilib- 
rium conditions for a circular plasma ring proceed- 
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Fig. 14. 

Cross section of a plasma ring. 



ing from the equations of magnetic hydrodynamics. 
In the first approximation the cross section of the 
plasma ring can be assumed to be the circle of ra- 
dius a. The plasma pressure and the current den- 
sity vanish outside the circle of radius a. The toroi- 
dal surface of the radius a is the boundary surface 
of the plasma ring which has no contact with the 
walls of the vacuum chamber in which it has been 
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generated (Fig. 14). This surface is also one of the 
magnetic surfaces of the toroidal trap. Another 
fundamental geometric parameter of the system is 
the major radius R of the plasma ring. Below we 
shall always assume that a is much less than R as is 
invariably the case in experimental environments. 
Note that when the ratio alR is small, the charac- 
teristics of the plasma ring are close to those of an 
infinite plasma cylinder. Since the problem has 
two fundamental geometric parameters a and R 
for two degress of freedom of the plasma ring, we 
have to find two equilibrium conditions — for the 
minor radius and for the major radius. For the 
first one we can neglect the effect of the toroidal 
structure for the minor radius. Thus, we can use 
the equilibrium equation for the plasma cylinder. 
Since the plasma in this case is under the effect of 
two fields, the equilibrium equations can be synthe- 
sized from Eqs. (61) and (62) which were derived 
for either the longitudinal field or the field of the 
current. We can derive such equilibrium equations 
for the plasma cylinder assuming that both the 
fields and the plasma pressure depend only on r, 
that is, on the distance between the given point 
and the axis of the plasma cylinder 

It can be easily shown that with the above con- 
dition, equilibrium necessitates the identity be- 
tween the sum of the plasma pressure p and the 
pressure B% (r)/8rt of the longitudinal magnetic 
field in the plasma averaged over the cross section 
of the plasma column and the sum of the pressures 
of the longitudinal field and the field of the current 
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at the boundary of the plasma column: 

p + Bl (r)/8n = Bl (o)/8n + B\ (a)/8n (64) 
Eq. (64) can also be written in the following form: 



2na?p = P/c* + [Bl (a) - B 2 e (r)]/4 (65) 

It can be demonstrated by detailed theoretical 
analysis that Eq. (65) is also valid when we take 
into consideration in the first approximation the 
toroidal structure of the system. The resulting 
equilibrium condition differs only by the corrected 
parameter B 9 (a), which, under rigorous treatment, 
is the longitudinal field at the surface of the plasma 
for (p = ji/2 (see Fig. 14). 

To derive the equilibrium condition for the 
major radius, we have to analyze the character and 
magnitude of forces which can result in variation 
of R. Let us restrict ourselves here to a rough 
qualitative analysis of these forces. There are 
three various factors giving rise to such forces. 

(1) Electrodynamic Stretching of the Ring Cur- 
rent. According to the fundamental laws of electro- 
dynamics, the ponderomotive force acting on a 
conductor carrying current always has such a direc- 
tion that it tends to increase its inductance. In our 
special case the inductance increases with the 
major radius R. 

Thus, we encounter here tensile stresses. This is 
almost self-evident, though. Since the pondero- 
motive force is due to the interaction between the 
current and its magnetic field it is proportional to 
/ 2 for given geometric parameters. 
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(2) Stretching of the Toroidal Ring under the 
Action of the Internal Pressure. To find the tensile 
force due to the above factor, we have to calculate 
the work performed by the plasma pressure for 
infinitesimal variation of R and to divide the 
result by 6R. The total tensile force proves to be 
proportional to a 2 p. 

(3) Radial Ponderomotive Force due to the Diffe- 
rence Between the Intensities of the Longitudinal 
Magnetic Field Inside and Outside the Plasma. When 
Bq (a) is higher than B% (a), this force is directed 
towards increasing R. 

The evaluation of the net effect of the above forces 
taking into account the equilibrium conditions for 
the minor radius shows that the net effect is to 
increase the radius R, that is, the pinch carrying 
current must stretch. This means the lack of equi- 
librium. But one step more and we reach our goal. 
The tensile force can be compensated by applying 
the magnetic field B± perpendicular to the equato- 
rial plane of the plasma ring (see Fig. 14). If we have 
appropriately chosen the direction of B±, it will 
act with the force 2nRIB±/c on the ring carrying 
the current / (referred to the total length of the 
ring). This force can compensate the net effect 
of the tensile forces discussed above. 

In the experimental machines in which toroidal 
plasma rings are generated, the tensile forces are 
compensated automatically, since the vacuum cham- 
ber has a thick conducting (metal) body. A dis- 
placement of the centre of the cross section of the 
plasma ring with respect to the centre of the cross 
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section of the toroidal conducting chamber body 
gives rise to eddy currents in the body. Eddy currents 
generate the transverse magnetic field which com- 
pensates the tensile forces. It can readily be shown 
that the compensating compressive force is pro- 
portional to I 2 and the displacement 6. Hence, 
equilibrium is reached at a certain value of 6. 

V. Shafranov has analyzed this equilibrium pro- 
blem theoretically and found the equilibrium dis- 
placement 6*: 

8 -£{"4+(<-£)ff^+VJ} « 

Here b is the radius of the cross section of the con- 
ducting body, B v (a) is the magnetic field generated 
by the current at the boundary of the plasma ring, 
and L t is the internal inductance of the plasma 
ring per unit length. Eq. (66) determines the posi- 
tion of the boundary magnetic surface of the plasma 
ring. Since the system is toroidal, the cross sections 
of the magnetic surfaces inside the plasma ring 
must not be represented by concentric circles (as is 
done for the plasma cylinder). Plasma carrying 
longitudinal current which is inside any of the 
internal magnetic surfaces has the same equilib- 
rium conditions with respect to the outer plasma 
layers as the conditions discussed above for the 
plasma ring surrounded by the conducting chamber 



* A rigorous definition of 6 is the distance 
between the geometric centre of the cross section 
of the plasma ring and the geometrical 
centre of the cross section of the body. 
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Fig. 15. 

Displacement of magnetic surfaces in a torus. 

body. Therefore, we expected the centres of the 
circles representing the cross sections of the magne- 
tic surfaces to be displaced outwards, towards the 
increasing R the larger, the smaller their radiuses 
(see Fig. 15). The largest displacement with respect 
to the centre of the cross section of the conducting 
body corresponds to the magnetic axis of the 
system where the current field B^ vanishes. 
These are the main results of theoretical analysis 
of equilibrium of the plasma rings carrying cur- 
rent. 

Let us note an interesting fact. The longitudinal 
magnetic field does not have any effect on Eq. (66) 
which determines the equilibrium position of the 
plasma ring in the chamber. It is as if this field 
is non-existent. It seems that logically the next 
step is to get rid of the longitudinal magnelic field 
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in the ring system altogether. This is quite admis- 
sible in the framework of the equilibrium theory. 
But once we become concerned with the stability 
of the plasma ring, the longitudinal field becomes 
again indispensable. 

When there is no longitudinal magnetic field, the 
bending and bottleneck-type deformations can de- 
velop without limitations so that in a very short 
time (of the order of a/v t , where v- t is the thermal 
velocity of ions) the regular shape of the plasma 
ring will spread out and the plasma will rapidly 
cool down owing to intensive interaction with the 
chamber walls. A strong longitudinal magnetic 
field produces a sort of a rigid frame of the field 
lines which stabilizes the plasma ring. Thus, even 
though the longitudinal field is indifferent with 
respect to the equilibrium conditions, it is instru- 
mental in making the plasma ring stable. This 
point, clearly, has to be discussed in more detail 
as is done below. 

Plasma is a good conductor; hence, when it tra- 
vels rapidly in a magnetic field and crosses its 
lines, induced currents appear which change the 
distribution of the intensity of the magnetic field. 
The general effect of this change is as if the plasma 
in its motion carries with it the field lines which 
behave as if they were "frozen" into the plasma or 
"glued" to it. To get a better understanding of the 
"freezing" of the field lines, let us consider the 
following special case of plasma motion. 

Let a magnetic field parallel to the axis z be 
applied to a volume of plasma. The field varies as 
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a given function in the plane xy perpendicular 
to the axis z. The plasma density is constant along 
the field lines but can vary with the coordinates x 
and y. Assume that the plasma travels perpendicu- 
lar to the held lines. Let us consider a small plasma 
column parallel to B whose cross section in the 
plane xy is dS. During the motion of plasma, this 
plasma column will be compressed or stretched so 
that its cross section will be changed. Since we 
deal here with rapid processes when diffusion 
can be neglected, the total number of particles 
in the column must remain constant. Hence, we 
obtain 

ndS = const (67) 

On the other hand, the magnetic flux through 
the cross section dS must remain constant, too, if 
the motion is rapid. If the plasma conductivity 
is high enough, even a slight variation of the mag- 
netic flux, if it occurs rapidly, gives rise to a high 
induced current. The magnetic field of this cur- 
rent compensates variation of the magnetic flux. 
Conservation of the flux means that 

BdS = const (68) 

Comparison of Eqs. (67) and (68) shows that under 
these circumstance the ratio Bin remains constant 
with the motion of plasma. If the plasma is com- 
pressed, the magnetic field in it increases and if 
plasma is expanded, the field decreases. Redistri- 
bution of the plasma density due to rapid deforma- 
tions results in appropriate variation of spatial 
distribution of the magnetic field. The "freezing" 
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of the field lines is just the effect of this relation- 
ship between deformations of plasma and varia- 
tions of the field which follows from conservation 
of the magnetic flux in a travelling volume ele- 
ment of the plasma. 

It should be noted, however, that the ratio Bin 
is not always constant; this is so only in the specific 
case discussed above. In different geometric condi- 
tions the relationship between the magnetic field 
and the density in the rapidly travelling plasma 
can differ from the above one. For instance, in an 
axially symmetric plasma jet travelling under the 
effect of electrodynamic forces due to current in the 
plasma, the parameter that remains constant for 
each plasma element during the motion is Blnr, 
where r is the distance between this element and 
the axis of the system. 

There is another noteworthy point. The "freezing" 
of field lines is not an effect which is unique to 
a plasma. A similar effect can be found in any 
good conductor travelling in a magnetic field with 
high enough speed. 

The stabilizing effect of the longitudinal field on 
the plasma ring is due to the "freezing" of field lines 
into the plasma. If a part of the ring gets more nar- 
row, the longitudinal field in this part increases 
owing to conservation of its flux; hence, the magne- 
tic pressure in the plasma in this part also increases. 

This force prevents the development of deforma- 
tion. Bending deformations result in elongation of 
the "frozen" lines of the longitudinal field and hence 
in increase of the Maxwellian stresses. This effect 
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can also prevent the development of deformations. 
The effect of the longitudinal magnetic field on 
some types of deformations of the plasma ring is 
a special case of those effects which are analyzed 
by the general theory of stability of plasma con- 
figurations. Later we shall discuss some of the 
general aspects of the theory but until that we 
shall sometimes note some results obtained by ap- 
plying the theory of stability to the specific problems 
we are interested in. 



THE THEORY OF STABILITY 

This theory can be applied, for instance, to 
selecting a particular technique for confinement and 
thermal insulation of the plasma ring using the 
combination of the longitudinal magnetic field 
and the field of the plasma current. As~predicted 
by the theory, the system has the highest stability 
with respect to the most dangerous large-scale 
deformations when the longitudinal magnetic field 
Bq is much higher than the current field B v . The 
stability condition can be written in the following 
form: 

(B e /B v ) (r/R) > 1 (69) 

Here r is the distance from a given point inside the 
pinch to the magnetic axis at which the field B v 
vanishes. The stability condition (69) has been 
derived independently jby M. D. Kruskal and 
V. D. Shafranov. Condition (69) must be satisfied 
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throughout the plasma ring. If we take into account 
the character of distribution of the current density 
over the cross section of the plasma ring (this distri- 
bution is always bell-shaped), we find that inequa- 
lity (69) is satisfied at any point of the plasma if it 
is satisfied with a margin at the plasma boundary, 
that is, for r = a. Hence, the Kruskal-Shafranov 
condition can be written in the following form: 

(B B /B V ) (a/R) > q (70) 

The quantity q is known as the margin of stability. 
It seems that stability is good enough when q is 
more than 3 unless the current distribution is very 
unusual. 

For a fairly long time, the experimenters in the 
USSR have been studying toroidal systems, where 
plasma is confined by the magnetic field of the 
current circulating in the plasma and the plasma 
ring is stabilized by a very high longitudinal 
field. Machines of this type are known as "Toka- 
maks". Fig. 16 shows schematically a tokamak-type 
installation. The toroidal discharge chamber encir- 
cles the core of the transformer which induces 
current in the plasma ring. The chamber has an 
inner shell made of thin-sheet stainless steel and 
an outer shell of thick-sheet copper. There is an 
electric insulation between two shells and the 
space between them is evacuated. Before each expe- 
riment, long-time heating is used to degas the 
inner shell, the so-called "liner", inside which the 
plasma is generated by the induced electric discharge. 
The outer copper shell stabilizes the plasma 
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Fig. 16. 

General view of a Tokamak. 



by eddy currents generated with displacement of 
the plasma ring. The longitudinal field B e is gene- 
rated by the magnetic coils installed on the cham- 
ber. In a medium-scale Tokamak (T-4), the radius 
a of the cross section of the plasma ring can be as 
high as 15-18 cm for R = 100 cm (the maximum 
value of a is determined by the size of the diaphragm 
inside the liner). 

The experiments with the T-4 Tokamak were 
carried out with the highest longitudinal field of 
about 37 kOe and the current density varied from 
100 to 200 A/cm 2 . The duration of the current 
pulse could be varied from a few milliseconds to 
50-60 ms. The mean plasma density in these expe- 
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riments was varied from n e « 1-10 13 to 4-10 13 cm -3 
by changing the initial gas pressure in the chamber. 
The experiments with other Tokamaks were carried 
out in a wider density range (from 2-10 12 to 
1-10 14 cm- 3 ). All the principal experiments were 
done with hydrogen and deuterium. 

The plasma in Tokamaks is heated naturally due 
to the Joule heating of the current flowing in it. 
Under these circumstances, the current density 
in the plasma is the principal factor determining 
the heating of plasma. The Kruskal-Shafranov sta- 
bility condition shows that the permissible current 
density increases as the longitudinal field. Hence, 
the higher B e the higher the plasma temperature 
th it can be obtained. 

The experimental results from Tokamaks fully 
support the validity of the Kruskal-Shafranov 
condition. The results show that if the stability 
margin q is sufficiently large (over 3 or 4) all the 
signs of the large-scale instabilities disappear and 
the plasma can be heated up to a high temperature. 
The current experiments with the T-4 machine 
with a current pulse of 25-60 ms and the longitudi- 
nal field of 37 kOe produce the hydrogen plasma 
with the following optimum parameters: n e « 
« 4-10 13 cm- 3 , Tgfal-iO 7 degrees, T, « 5-10 6 
degrees. Analysis of the energy balance of the plas- 
ma yields energy liberated in the plasma ring owing 
to the Joule heating and the thermal energy flux 
from the plasma ring. The magnitude of this flux 
is a gauge of the efficiency of thermal insulation of 
plasma which is related to the magnetic field. The 
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efficiency of thermal insulation can be characte- 
rized by the parameter known as the mean time 
of energy retention in the plasma and is defined 
in the following way 

t e = WIQ (71) 

Here W is the total energy stored in the plasma 
ring, Q is the energy leaving the plasma per unit 
time (lost power) and T B is the mean time of energy 
retention. The time t E in Tokamaks can be as 
long as 45 ms under the optimum conditions (high 
longitudinal field, high plasma current, and high 
enough plasma density). 

The mean lifetimes of fast electrons and ions in 
the plasma can also be determined in experiments. 
In the T-4 machine, these lifetimes are more than 
20-30 ms under the optimum conditions. An ion 
travels about 10 km in the chamber during this 
period and undergoes about a hundred Coulomb 
collisions. A comparison of the confinement time 
for particles in Tokamaks and the time predicted 
from the diffusion theory using the expression for 
the diffusion coefficient in the magnetic field B e 
given above shows that the experimental confine- 
ment time is by a few orders of magnitude smaller 
than the predicted time. A few years ago this fact 
was assumed to indicate the existence of "anoma- 
lous" diffusion due to development of a "slow" 
instability of plasma. This idea was supported 
by the measured heat losses from the plasma ring. 
These losses were by a few orders of magnitude 
higher than those predicted by the theory of plasma 
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heat conductivity in a magnetic field using in the 
expression for the heat conductivity coefficient 
the value of B equal to the longitudinal field in 
a Tokamak. 

It would seem that the anomalous heat conduc- 
tivity can be attributed to the development of 
"slow" instabilities. We cannot be quite sure, 
however, that high loss rates of plasma energy 
and particles in Tokamaks are due to any insta- 
bilities. This is because rigorous classical theory 
of diffusion and heat conductivity shows that the 
diffusion and heat conductivity coefficients in 
Tokamaks should be higher by about three orders 
of magnitude than it was assumed earlier. Intui- 
tively, it can readily be seen from {"qualitative 
considerations that the formulas for coefficients 
of diffusion and heat conductivity for the Tokamak 
magnetic systems should be modified. Particles 
in Tokamaks travel in a composite magnetic field 
and the longitudinal field B e alone cannot hold 
them. Hence, the current field B^ must play the 
major part in the transport phenomena and, since 
By <^ B e , the transport rates should be expected 
to be higher than in the uniform field B e . 

The correct quantitative theory of transport 
phenomena in toroidal systems was developed in 
two steps. At first, B. B. Kadomtzev and V. D. Shaf- 
ranov took into consideration the effect of toroidal 
drift on the trajectories of particles and found 
that the coefficients of diffusion and heat condu- 
ctivity must be increased by a factor of about 2g 2 , 
where q = \BJB^ (a)] (a/B) if the plasma density 
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Fig. 17. 

Trajectory of a confined particle. 



is high enough and there is no physical meaning in 
classifying particles as through-going and confined. 
Thus, the diffusion and heat conductivity coeffi- 
cients under normal conditions of the experiment 
must be higher by a factor of 20-30. 

The next theoretical advance was done by R. Z. 
Sagdeev and A. A. Galeev. They noticed that a 
toroidal drift had an especially strong effect on 
the confined particles which oscillated along the 
lines of force being reflected from the high-field 
regions. Fig. 17 shows such a trajectory and its 
projection on the plane in which the system's axis 
lies and which rotates by the angle with the particle. 
The particle also oscillates along the radial direction 
so that the trajectories of particles in the magnetic 



91 



field are mixed and the rates of diffusion and heat 
transfer increase greatly. These rates are higher 
than the earlier predictions (B. B. Kadomtzev 
and V. D. Shafranov) by a factor of (R/a) 3 ' 2 , that 
is, about 30 times under the experimental condi- 
tions. This effect occurs when collisions in the 
plasma are rare enough; therefore, we can clas- 
sify the particles into confined and through-going 
ones. In Tokamak experiments, collisions are 
typically rare, hence, the coefficients of diffusion 
and heat conductivity must be higher by three 
orders of magnitude than assumed earlier. 

It is yet early to assert that the experimental 
results are in a very good agreement with the modi- 
fied classical theory. The current experimental and 
theoretical accuracies of the values of the diffu- 
sion and heat conductivity coefficients are such 
that if they differ by a factor of two or three, this 
does not necessarily mean contradiction in the 
results. 

The available results suggest that at least under 
some experimental conditions when a plasma has 
a relatively high ion temperature and a not-too-low 
density, the measured rates of diffusion and heat 
transfer are close to those predicted by the theory 
of R. Z. Sagdeev and A. A. Galeev. 

However, we still cannot rule out the possibility 
that the slowly developing so-called "drift" insta- 
bilities also contribute to the loss of particles and 
heat energy from a high-temperature plasma. But 
even now it is clear that the plasma temperature 
and the time of energy retention increase conside- 
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rably with further increase in the magnetic held 
and the cross section radius of the plasma ring in 
Tokamaks. 

In experiments with toroidal systems, the pri- 
mary aim remains to be the same — to increase 
the plasma parameters (temperature, pressure and 
lifetimes of particles) further. Only when stable 
high-temperature plasma configurations have been 
generated, we shall be able to study in earnest 
the fine physical effects typical of plasma, in par- 
ticular, those sophisticated tunes that could be 
readily played with unique resonator. 

10. 

PLASMA BEHAVIOUR " 

IN TRAPS WITH MAGNETIC MIRRORS 

Plasma is confined in such systems owing to 
reflection of particles from high-field regions. In the 
magnetic mirror traps the distribution of velocities 
of plasma particles is non-isotropic. 

Let the maximum intensity of the magnetic 
field along one of the field lines is 5 max . This value 
can be assumed to be the same at both ends of the 
trap since, if the values were different at the ends, 
the confinement of particles would be determined 
only by the reflection of particles from the mirror 
with a lower field. 

Let us consider a small segment of the field line 
at which the field intensity is B. For a long lifetime 
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Fig. 18, 

Loss cones in a trap with magnetic mirrors. 

of the plasma in the trap, only those charged par- 
ticles, for which the angle a between the velocity 
and the fie ld line is larger than arc sin a = 
= VB'Bmaxi can be within this segment. All 
the particles with smaller angles a will leave the 
trap within a very short period (of the order of 
llv t , where I is the trap's length and v t is the thermal 
velocity of ions). In particular, near the magnetic 
mirrors, where B is close to 5 max only such particles 
can be found whose longitudinal velocity is rela- 
tively small and the angle a is close to 90°. Fig. 18 
presents the sphere of all possible directions of the 
particle's velocity. The velocity vectors of the 
particles which freely leave the trap lie within 
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two cones in this sphere. These cones are called 
loss cones. In such a trap the plasma pressure, clear- 
ly, is not a scalar quantity. We have to distinguish 
between the longitudinal pressure 

Pl|r= niiriivli + m^n,) 
and the transverse pressure 

p L = n (rriivii + mj)\ e )i2 

If the velocity vector of a particle gets into a loss 
cone owing to a Coulomb collision, this particle 
will immediately leave the plasma. If the angle 
of the loss cone is large enough, practically one 
Coulomb collision is sufficient to expel a confined 
particle from the trap with magnetic mirrors. There- 
fore, the mean lifetime of particles in such traps 
must be fairly close by the order of magnitude to 
the mean time x ti between two successive ion-ion 
collisions. The ion-ion collisions play here a deci- 
sive role, since electrons, the collisions between 
which are more frequent, cannot leave the trap 
independently of ions. At the moment of plasma 
generation in the trap, a relatively small number 
of electrons will leave it so that the plasma will 
have a positive charge and the flows of particles 
of both signs will be balanced. 

The above results show that even under optimum 
conditions when the loss of particles is due only to 
Coulomb collisions the open-ended traps with 
magnetic mirrors have a much lower confinement 
times than the toroidal traps in which the mea- 
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sured confinement times are higher than T (j by two 
orders of magnitude. But an important experimen- 
tal advantage of the open-ended traps is the possi- 
bility of using a variety of techniques for generat- 
ing a high-temperature plasma. In open-ended traps, 
a plasma can be generated by injection of fast 
particle flows, capture of plasma jets, high-frequency 
heating of cold plasma, production of fast ions in 
cold plasma jet by specially generating instabili- 
ties, etc. These traps also make it possible to vary 
in wide ranges the energy spectra of electrons and 
ions and to generate a plasma with high-energy 
ions and low-energy electrons or, on the contrary, 
with very fast electrons and relatively cold ions. 
At an early stage of the experimental open-ended 
system studies of plasma, the only systems used 
were those with two magnetic mirrors at the ends. 
The magnetic field in such traps increases along 
the field lines in both directions from the central 
region where the plasma should be. At the same 
time, the field decreases radially. Numerous expe- 
riments have shown that this simplest open-ended 
trap is not suitable for long confinement of plasma 
with a hot ion component. For the plasma density 
over 10M0 8 cm 4 , the lifetime of plasma in a sim- 
ple two-mirror trap is not more than a few tens of 
microseconds. The plasma flows out from the 
trap owing to the development of the so-called 
flute instabilities. This instability is due to the 
fact that the plasma is a diamagnetic and therefore 
it can easily travel towards a weaker field by pro- 
ducing "tongues" which seep between the field 
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Fig. 19- 

Diagram of a magnetic system 
with a combined field. 



lines (between the "tongues" there appear flutes 
without plasma). 

To prevent the development of flute instability, 
in the simplest magnetic mirror traps such tech- 
niques have to be applied which by themselves 
result in considerable deterioration of thermal 
insulation of plasma which leads to its cooling 
(low vacuum, intense low-voltage arc in the space 
filled with plasma, etc.). 

In 1961 M. Ioffe was the first to experiment with 
an open-ended trap with the magnetic field increas- 
ing in all directions from the plasma region. Fig. 19 
illustrates the magnetic system of such a trap. The 
magnetic field is produced by coaxial coils and 
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Fig. 20. 

Field of linear conductors. 



six linear conductors carrying current which are 
placed symmetrically with respect to the central 
axis of coils. Currents flowing in the neighbouring 
conductors have.. opposite directions. This system 
without linear conductors would be a simple mir- 
ror trap with a radially decreasing field. The linear 
conductors generate the field whose geometry is 
shown in Fig. 20. This field varies as r 2 , where r 
is the distance to the symmetry axis. In combined 
fields of the coils and linear conductors, the result- 
ing field intensity increases in all directions from 
the central region of the trap where a plasma is 
generated. This is known as the minimum-S con- 
cept. 
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The experiments with minimum-i? machines have 
shown that they make possible stable confinement 
of plasma with a density up to about 10 u cm -3 
and the ion energy of the order of 1 keV (T t ~ 
— 1 • 10 7 degrees). The lifetimes of plasma particles 
in such machines amount to tens of milliseconds. 
Thus, a technique has been found for suppressing one 
of the most dangerous rapidly developing insta- 
bilities typical of hot plasma. Therefore, all modern 
open-ended hot-plasma systems are designed pro- 
ceeding from the minimum-/? concept. 

It is yet unclear whether higher temperatures 
and plasma densities in the minimum-5 machines 
can give rise to other types of plasma instability. 
Theoretical predictions give some food for these 
thoughts. Analysis of plasma behaviour in open- 
ended machines with magnetic mirrors shows that 
the peculiar distribution of particle velocites in 
such systems gives rise to various mechanisms 
which generate so-called kinetic instabilities. The 
kinetic instabilities are due to amplification of the 
plasma oscillations by directed flows of ions or 
electrons. These instabilities can readily develop 
if the velocity distribution function has a maximum 
at a nonzero velocity. The existence of such a maxi- 
mum is equivalent to the existence of a directed flow 
of particles. 

If the plasma particles have the Maxwell velocity 
distribution function (the dotted line in Fig. 21), 
the oscillations are not amplified. But in an open- 
ended machine, the velocity distribution function 
for confined particles must have a distinct maximum 
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Fig. 21. 

Velocity distributions of particles: 
Maxwell distribution— dashed line; 
distribution in the presence 
of a loss cone — solid line. 



for the transverse velocity component owing to 
the existence of the loss cone (see the solid line 
in Fig. 21). With this shape of velocity distribu- 
tion, there can occur amplification of ion-electron 
oscillations resulting in a sharp increase of the 
plasma loss from the trap. 

The first to analyze theoretically this type of 
instability were Rosenbluth and Post; they have 
shown that this instability develops when the 
length of the plasma clot trapped in the machine 
is more than a few hundreds Larmor radiuses. 
This instability can, apparently, be stabilized 
in a short magnetic system. 
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The theory predicts another type of instability 
which is due to non-uniform distribution of plasma 
density in the direction perpendicular to the field 
(this non-uniformity is always present merely because 
of the finite cross section of the plasma configu- 
ration). The drift motion of ions in non-uniform 
plasma can give rise to amplification of ion oscilla- 
tions of a frequency close to the Larmor frequency 
co Bj . The oscillations are amplified if the Larmor 
radius of ions amounts to about a lB 2 /(4nnm i c 2 )] 2/3 , 
where a is the radius of the plasma. 

The experimental plasmas in open-ended ma- 
chines have the parameters close to values at which 
the conditions for the development of the above 
instabilities are satisfied. However, up till now 
such instabilities have not been observed. There- 
fore, we can hope that we can attain higher tem- 
peratures and densities for stable plasma confined 
in the minimum-5 open-ended machines. 

11. 

STABILITY OF PLASMA 
CONFIGURATIONS 

Even a very brief discussion of plasma physics 
cannot ignore the central problem of this science, 
namely, the stability of plasma configurations. 
The advances in stability studies determine the 
prospects for increasing the temperature and life- 
time of hot plasma. 
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Our general discussion of plasma stability will 
lack mathematical rigour. With this in mind, let 
us turn to our subject. Let us assume that a con- 
fined plasma configuration is suspended in vacuum 
by a magnetic field and is at magnetohydrodyna- 
mic equilibrium. In such a plasma, the temperature 
and density remain constant and Eq. (52) is satis- 
fied at each point. Will the equilibrium state be 
preserved for a long time or the random fluctuations 
in the plasma will generate increasing instabilities 
and the plasma will rapidly spread out over the 
available volume? This is the generally formulated 
stability problem. A significant factor in this 
problem is that different instability mechanisms 
produce perturbations with widely different rates 
of increase in the plasma. 

The development of theory should, to a certain 
extent, take this fact into account. Thus, the 
origin of instabilities and the techniques for their 
stabilization should be considered, starting from 
the instabilities which give rise to rapid macrosco- 
pic motion of plasma and are, hence, the most dan- 
gerous. When these, the most significant instabili- 
ties have been analyzed, we can then analyze less 
destructive instabilities and so on. It can appear 
that if we take into consideration all the types of 
instability of plasma confined in a magnetic field, 
we shall find that an absolutely stable plasma con- 
figuration is unfeasible. But absolute stability is, 
maybe, not essential. In practice, we must learn 
how to inhibit the development of various insta- 
bilities, so that non-stabilized weak and slow-deve- 
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loping instabilities would not seriously hinder 
accumulation and long-time retention of thermal 
energy in high-temperature plasma. 

The most dangerous instabilities are those which 
give rise to motion of macroscopic plasma elements 
with the velocities of the order of thermal ion velo- 
city. For such a rapid motion, the plasma in the 
magnetic field can be fully described by the model 
of an ideally conducting liquid. Therefore, rapidly 
developing large-scale instabilities are known as 
the magnetohydrodynamic instabilities. They de- 
velop owing to the diamagnetic properties of plas- 
ma. The plasma tends to shift towards a decreasing 
field; therefore, if the plasma is in a region where 
the field decreases in the outward direction from 
the plasma boundary, this boundary can prove to 
be unstable. As discussed above, in open-ended 
macb'nes, we can generate the minimum-5 mag- 
netic field which provides for the magnetohydrody- 
namic stability of the plasma in the machine. 

For the closed traps, however, it has been estab- 
lished that it is impossible to generate a magnetic 
field increasing outwards from the plasma boun- 
dary at each surface point of a toroidal plasma 
configuration. The grad | B | component normal 
to the plasma surface has different signs at diffe- 
rent points of the boundary. For instance, in Toka- 
maks the value of B decreases outwards at the 
outside of the plasma ring and increases at the insi- 
de of the ring. Thus, a question arises whether 
the plasma can flow towards a decreasing field in 
the form of individual "tongues". To answer this 
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Fig. 22. 

"Tongue" on the plasma surface. 



question, we have to differentiate between a high- 
pressure plasma for which p ~ B 2 /8n and a low- 
pressure plasma for which p <^ B 2 /8n. On the 
surface of the high-pressure plasma, there can 
develop local perturbations in the form of "ton- 
gues" (Fig. 22). Since the field is "frozen" into the 
plasma, the development of the tongue results in 
an increase of magnetic energy and in curving of 
field lines. An appropriate work will be done by 
the expanding high-pressure plasma at the expense 
of its thermal energy. If the tongue extends towards 
a decreasing field, it will extend further so that 
the plasma boundary will be unstable. In this case 
instability has a local nature since it is determined 
by the local field geometry. In toroidal traps there 
always exist parts of the plasma ring where the 
field geometry encourages the development of 
tongues. Hence, the high-pressure plasma is un- 
stable in such systems. 

If p <^ S 3 /8jt as is the case in modern toroidal 
machines (for instance, Tokamaksand Stellarators), 
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Fig. 23. 

Development of flutes on the plasma surface. 

plasma perturbations cannot distort noticeably 
the shape of the field lines. The plasma just lacks 
an extra thermal energy to do this. Hence, small 
local perturbations of tho tongue-type are stabili- 
zed automatically and the only! perturbations inside 
the plasma or at its boundary that can occur are the 
rearrangements of the whole field tubes and the 
formation of flutes (Fig. 23). 

The plasma which filled the field tube formed 
by a thin bunch of field lines tends to expand and 
will, therefore, migrate towards the side where 
the tube volume will increase. This volume is 

\ 6Sdl, where 8S is the area of the"' tube cross section 

and dl is the elementary length of the field line. 
Since the magnetic flux 6d> is constant along the 
length of the tube, we have 

8v = j 6SBdl/B = j SOdl/S = 80> J dl/B (72) 

The flux 8<t> remains constant for all movements of 
the tube. Hence, the tube" volume varies as \ dl/B. 
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Since a plasma, as any other gas, tends to increase 
its volume in motion of a tube with the plasma, 

the parameter u = — \ dllB plays a part similar 

to that of the potential energy. Movements of 
individual plasma elements with the exchange of 
the field tubes are known as rearrangement or con- 
vective deformations. When they reach the plasma 
boundary with an external field, the plasma surface 
acquires a fluted structure oriented along the field 
lines. These are the "flute" instabilities' briefly men- 
tioned above. 

Making use of the concept of "potential energy" 
u (this is a purely arbitrary term), we can formu- 
late for'the low-pressure plasma thestability criterion 
with respect to flute instabilities. The plasma boun- 
dary is stable if u increases along the field tube with 
the plasma towards its intersection with the plasma 
boundary, that is, if the following inequality is 
satisfied: 

6 \ dllB < (73) 

The variation of the integral is taken along the 
perpendicular to the plasma surface between two 
infinitesimally close field lines. If the inequality 
is not satisfied, the boundary is unstable. 
r3The stability condition (73) means physically 
that the necessary and sufficient condition of stabi- 
lity is the increase of the averaged field along the 
field line outwards from the plasma boundary. 
This is a modified minimum-5 concept. Note that 



106 



according to Eq. (72) [ dl/B -■= dv/b<t>. Therefore, 

the quantity \ dl/B can also be named the specific 

volume of the field tube. 

To be able to use the above stability criterion, 
we have, at first, to eliminate the uncertainty in 
the expression for u, where no limits of integration 
along the field line are given. There is no uncer- 
tainty when the field lines are closed at the plasma 

surface. In this special case, the integral \ dl/B, 

clearly, must be taken along the length of the 
field line. We can introduce here the following 
improved definition of the potential energy for the 
magnetic surface with closed field lines: 

u = -UN <J> dl/B (74) 

Here N is the number of turns around the magnetic 
system after which the field line is closed upon 
itself. This can be naturally extended to the case 
of closed field lines: 

u = lim UN \ dl/B (75) 

Here N and the length of integration tend to infi- 
nity. The parameter u can be shown to be a defi- 
nite characteristic of the magnetic surface. There- 
fore, the expression for u can be rewritten to include 
only the parameters of the given magnetic surface 
so that any uncertainty in the length of integration 
is eliminated. 
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The volume V corresponding to a closed magnetic 
surface can be treated as a function of the flux tp e 
of the longitudinal magnetic field passing through 
its cross section: V = F(<J> e ). Using this notation 
we can write the potential energy as u —- —V (<J>e), 
so that the stability condition can be written in 
the following form: 

F"(<D 9 )<0 (76) 

The form (76) of the stability condition is valid 
not only for the systems in which the field lines 
are closed, but also for the general case when the 
field lines on the magnetic surface are infinite. 

For each plasma configuration we can plot a 
"magnetic well", that is, the curve u (p), where p 
is the distance between a given surface and the 
magnetic axis (at a small distance from the axis 
we have O e ~ p 2 ). The relative depth of the mag- 
netic well is 

Au/|u | = [»(p ) -a(0)]/|u(0) |! 

where u (p ) and u (0) are the values of u at the 
plasma surface and at the magnetic axis. The 
higher the parameter Au/\ u | the better the sta- 
bility of plasma with respect to the convective 
instabilities. If this parameter is negative, that is, 
the system has a magnetic "peak" instead of the 
well, the plasma is unstable. Analysis of the fun- 
ction u (p) (or the equivalent function V (<J> e )) 
shows that the existence or lack of the magnetic 
well is determined by the arrangement of magnetic 
surfaces with respect to the magnetic axis. 
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Fig. 24. 

Displacement of magnetic surfaces inwards (a) 

and outwards (6) 

from the magnetic axis. 

Fig. 24a and b illustrates two different arrange- 
ments of magnetic surfaces in toroidal machines: 
at the left diagram the cross sections of magnetic 
surfaces are shifted inwards from the magnetic axis, 
that is, towards the geometrical centre of the 
system as a whole and the displacement is the 
larger the larger is p; at the second diagram the 
situation is reversed and the magnetic surfaces are 
displaced outwards from the magnetic axis. In the 
first case u increases with p; hence, there is a mag- 
netic well and the plasma is stable. In the second 
case there is no magnetic well and the plasma is 
unstable. These results follow from the rigorous 
theoretical analysis but qualitatively they are 
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almost self-evident. If the magnetic surface is 
displaced inwards, it gets into a stronger longitudi- 
nal field. Therefore, the flux increases faster than 
the volume and V" (O e ) must be negative, which 
means that the plasma is stable. When the magne- 
tic surfaces are displaced outwards, the effect 
must be reversed. 

Let us illustrate the above discussion by the 
examples of Tokamaks and Stellarators. As mentio- 
ned above, in Tokamaks the magnetic surfaces 
are displaced inwards with respect to the magnetic 
axis. Therefore, the plasma ring in Tokamaks must 
be stable with respect to convective magnetohydro- 
dynamic instabilities and, in particular, to flute 
instabilities at the surface. However, in the existing 
Tokamaks the relative depth of the magnetic well 
is very small (of the order of a few per cent). 

Analysis of the geometry of magnetic surfaces 
in stellarators shows that magnetic surfaces in 
them are displaced outwards. Hence, we have 
V" (^e) > an d the plasma ring can be unstable. 

12. 

SHEAR STABILIZATION 

The magnetohydrodynamic instabilities can be 
suppressed not only by means of a magnetic well 
but also by the so-called shear, that is, crossing 
of the field lines. If there is the shear in a system 
where the field line rotates around the magnetic 
axis, this means that the rotation angle of the field 
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lines is a function of r, that is, it varies from the 
magnetic axis to the periphery of the plasma ring. 
The rotation angle depends on the geometry of the 
magnetic system in the presence of a plasma ring. 
For instance, in Tokamaks it can readily be seen that 
for alR <<^ 1 the angle of rotation of the field line 
around the magnetic axis is (per unit length) 

u- = BJrB e 

The shear exists when \i is a function of r, that is, 
if d\xldr is non-zero. To compare shears in diffe- 
rent plasma systems we have to make use of the 
dimensionless parameters. In the publications on 
the stability theory the usual shear parameter is 

8, = r 2 d(i/dr (77) 

In the existing experimental machines the shear 
B s is not more than a few per cent. 

It is widely believed that the shear is some uni- 
versal technique for suppressing a very extensive 
range of plasma instabilities. Let us attempt illu- 
strating the shear stabilization by discussing the 
magnetohydrodynamic instabilities. 

Assume that we have a plasma with a low pressure 
8np/B 2 <^ 1 . If the shear 8 is close to zero and 
there is no magnetic well, then the convective 
instabilities can freely travel along the radius by 
rearrangement of the field tubes with the plasma 
and come out on the surface of the plasma ring 
in the form of flute instabilities. With an increas- 
ing shear such instabilities will be stabilized, since 
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a radial displacement of plasma in this case will 
result in a strong distortion of the field. When 
9 S =^0, then a deformation which had at a certain 
surface the form of a protrusion coinciding with 
the field line must curve at the neighbouring magne- 
tic surface since the field lines have rotated. The 
result must be a plasma tongue of a complicated 
shape. Its edges are as though fixed along the inter- 
nal field line and the radially farthest ridge is paral- 
lel to the external field line. The extra magnetic 
energy required for this distortion of the field's shape 
can be obtained only at the expense of the work of 
plasma pressure (remember that the field lines are 
frozen into the plasma). Hence, with a finite shear, 
perturbations in the plasma can be stabilized if 
the ratio 8np/B 2 is small enough. It can happen, 
however, that with a low shear a hot plasma can 
be stabilized only at a very low density. 



13. 

OTHER PLASMA INSTABILITIES 



The above discussion dealt with the rapidly 
developing large-scale instabilities for which plas- 
ma is an ideal conductor. The real plasma has a 
finite conductivity. Therefore, for slowly develop- 
ing instabilities, the field lines cannot be assumed 
to be frozen into the plasma and are not carried 
with it. The "thawing out" of the field lines must 
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facilitate propagation of perturbations in the plas- 
ma and can make the effectiveness of the above sta- 
bilization mechanisms insufficient. This extends 
the range of possible plasma instabilities. This can 
be illustrated by the so-called current-convective, 
or "corkscrew" instability which develops in the 
plasma when there is a current parallel to the mag- 
netic field. This type of instability has the follow- 
ing origin. A temperature fluctuation in a plasma 
with a longitudinal current changes the electric 
field in the plasma owing i.o the temperature de- 
pendence of conductivity. Under certain condi- 
tions, this can give rise to a drift motion perpendi- 
cular to B which involves the entire plasma layers. 
If the plasma has a temperature gradient, the drifts 
can, in their turn, amplify the initial thermal 
fluctuations owing to the replacement of colder 
layers by hotter ones. 

The instabilities due to finite conductivity of 
the plasma are usually known as dissipative insta- 
bilities. The perturbations with these instabilities 
are large in scale but have relatively low rise incre- 
ments. The dissipative instabilities can decrease 
the efficiency of plasma heating owing to increas- 
ing heat losses. But with increasing plasma tempe- 
rature, their significance decreases (since the plasma 
conductivity increases with temperature). 

The effect of dissipative instabilities on energy 
retention in a hot plasma seems to be practically 
confined to the electron temperature range below 
10 s degrees. When we have passed the initial heat- 
ing stage and generated a plasma with a higher 
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electron temperature, the dissipation processes 
practically cease to affect the plasma stability. 

If we have managed to produce a plasma confi- 
guration in which there is no development of mag- 
netohydrodynamic instabilities including the dis- 
sipative ones, it still does not follow that diffusion 
and heat conductivity in this plasma are due only 
to pair collisions of charged particles. There is 
another type of instabilities which can make the 
confined plasma "chronically sick" and result in 
its. premature "demise" following spreading out 
or cooling. These are the so-called drift instabili- 
ties. 

Before discussing the mechanism of drift insta- 
bilities, I would like to mention a problem of a more 
general nature. In a plasma there can develop and 
propagate a variety of oscillations and waves 
in a very wide frequency range. For instance, 
when electrons are displaced with respect to ions, 
this gives rise to high-frequency electron oscillations 
with an angular frequency (x) = y Anne 2 /m e . These 
oscillations can propagate in the plasma as longi- 
tudinal electrostatic waves. There can also be genera- 
ted waves of a much lower frequency due to the longi- 
tudinal oscillations of ions— the so-called ion 
sound. Along the field lines of the magnetic field, 
in the plasma there can propagate transverse ele- 
ctromagnetic waves (the so-called Alfven waves) 
whose velocity is 5/|/4jxn/Wj and is much lower 
than the light velocity when the plasma density 
is high. 
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Fig. 25. 

Interaction between a particle and a wave. 

We shall not go here into the details of wave 
problems but mention only one general question, 
namely, the interaction between waves and particle 
flows in the plasma. The character of this interac- 
tion will determine amplification or attenuation of 
waves of a given type. Assume that a longitudinal 
wave parallel to the axis x has been generated in 
the plasma. The wave process here is the periodic 
oscillations of the electric potential travelling 
with the velocity v t = co//e, where o> is the cyclic 
frequency and k is the wave vector (k = 2n/X). 

Let us analyze the interaction between the plasma 
particles and this running periodic potential bar- 
rier. Electrons or ions whose velocity component 
along the axis x greatly differs from v } will have small 
oscillations of velocity in the region of variable 
potential with the conservation of the mean energy. 
But if a particle has the velocity component v x close 
enough to v f , the particle and the wave will exchange 
energy. To understand this energy exchange, let us 
make use of the reference system which travels with 
the wave. Fig. 25 shows variation of the potential 
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energy of the particle in this system. Let the ampli- 
tude of oscillations of the potential energy be u m . 
The component of the kinetic energy of the particle 
along the axis x in the travelling reference system 
is m (v x — (&lk) z l2. If this energy is less than u m , 
the particle will be reflected by the potential bar- 
rier. We have to distinguish between two cases 
here. In the first case we have 

w/k + Y2uJm >v x > (a/k (78) 

and in the second case we have 

a/k > v x > <a/k — V 2u m /m (79) 

In the hrsl, case the particle 1 catches up with 
the wave and in the second case the wave catches 
up with the particle 2. Following the collision with 
the potential barrier, the velocity component v x of 
the first particle will reverse its sign and will be 
atlk — v x . In the laboratory reference system it 
will be 2co//« — v x . Eq. (78) shows that this velocity 
is less than v x . Thus, the particle's velocity has 
decreased following the collision, so that a part 
of its kinetic energy has been transferred to the 
wave. On the contrary, the collision of the second 
particle with the potential barrier will increase 
its velocity owing to the energy exchange with 
the wave. This simple analysis shows that the 
amplitude of the longitudinal electric wave in the 
plasma will increase or decrease depending on the 
majority of the plasma particles catching up with 
the wave or falling behind it. In other words, ampli- 
fication or attenuation of the wave will be determi- 
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Fig. 26. 

Velocity distribution of particles: 

a — Maxwell distribution; 

b — distribution in the presence of a beam. 



ned by the shape of the velocity distribution of 
particles for the velocity component parallel to 
the wave propagation direction. 

Fig. 20 presents two different velocity distribu- 
tions for a given direction (the first distribution 
corresponds to the Maxwell velocity distribution). 
Here we have df/dv x < in the whole velocity 
range. Hence, for any phase velocity of the wave 
and any amplitude of potential oscillations, there 
will be more particles that fall behind than those 
that catch up with the wave. Therefore, the wave 
will be attenuated. This mechanism is known as 
the Landau damping. The second velocity distri- 
bution function has the maximum at v x = v . To 
the left of this maximum, we have df/dv x > O.Hence, 
the waves with the phase velocities in this range 
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will be amplified owing to the energy transfer 
from the particles. 

Note that if the velocity distribution for the 
velocity component v x has a maximum at v x =^0, 
this can be interpreted as a result of a directed 
flow of particles in the plasma. (This means that 
the distribution in Fig. 266 can be represented as 
a superposition of a steadily decreasing function 
and an individual peak.) 

It is clear that the amplification of waves in 
the plasma is due just to the energy transfer from 
the particles of the directed flow. The energy trans- 
fer from a directed flow of electrons or ions to waves 
in the plasma decelerates the flow and increases the 
periodic electric field in the plasma. Various speci- 
fic effects due to this physical mechanism are known 
as the beam instabilities. They are encountered 
in various phenomena occurring in plasma. Howe- 
ver, the development of such an instability does 
not mean by itself that the thermal insulation of 
plasma has been disrupted and the diffusion rate 
has increased. The drifts contributing to the loss 
of particles and energy from the plasma can deve- 
lop only if the electric field of the wave due to the 
beam instability has a component perpendicular to B. 

This is just the case for the beam instabilities 
discussed below. These instabilities are, in a sense, 
universal, since they are due to the fact that the 
plasma has non-uniform density and temperature 
in the direction perpendicular to the magnetic field. 
Any confined plasma configuration'' is not uniform 
in this respect for the simple reason that near the 
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plasma boundary the concentration should sharply 
decrease and the temperature should also drop. 
Let. us consider the simplest plasma configura- 
tion—a long cylindrical column of a low-pressure 
plasma in the longitudinal magnetic field. We shall 
assume the magnetic field to be uniform. The equi- 
librium equation yields 

U = clB | grad p | (80) 

where j L is the density of diamagnetic current in 
the plasma. In the case of cylindrical symmetry, 
Eq. (80) can be written as 

1 1 (81) 

Here the current will be perpendicular to the pres- 
sure gradient and the vector B, that is, it will flow 
in the azimuthal direction. Since the diamagnetic 
current is due to the Larmor motion of electrons 
and ions, we can write 

7± =« e (u* — u r ) 

where u t and u e are the averaged drift velocities of 
electrons and ions. The drift velocities have also 
the azimuthal direction. We have here 



7x = clB 



Ui 



n a B 



dpi 



dr 



c I dp e 

Ue = 1~b ~~dT 



(82) 



If electrons and ions have close temperatures, so 
that we can introduce the general plasma tempera- 
ture T, we can write 

-' ■-■■' - ~ ~ ekT 

Xti r^s Ua r**/ 77— 

1 e eB 



p dr I eats v 
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Here a is the characteristic linear size of pressure 
non-uniformity in the plasma: 

p dr 

If the concentration drops from the axis to the 
boundary of the plasma column more sharply than 
does the temperature, we have 

_j 1 dn 

~" n dr 

If the temperature drop is more sharp, we obtain 

_. 1 dT 

a ^-T-dF 

Analysis results suggest that when in the plasma 
there develop periodic perturbations of the wave 
character which propagate at a small angle to the 
azimuth, the phase velocity of this wave to/ky is 
ckTIeaB and such waves are known as the drift 
waves. Since the drift azimuthal velocities of 
particles are close to the phase velocity of waves 
propagating in the same direction, the particles 
should be expected to interact strongly with each 
other so that conditions can be established for 
rapid amplification of oscillations. Increase in the 
amplitude of oscillations results in development 
of strong azimuthal electric fields, therefore, when 
the plasma has no uniform density and temperature, 
there must develop radial flows of heat and particles 
which shorten the lifetime of a confined plasma 
configuration. That is the effect of drift instabili- 
ties. 
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Detailed theoretical analysis shows that the ampli- 
fication of drift waves due to non-uniform plasma 
density is effectively inhibited by a high enough 
shear. The drift instability due to the temperature 
gradient is highly dangerous. Instabilities of this 
type are known as drift temperature instabilities. 
Such instabilities develop when d (\nT)/d (In n) > 1 
and considerably increase the heat losses from the 
plasma. Though an additional heat loss due to 
the drift temperature instability in closed toroidal 
systems can be sharply reduced by shear, it still 
can be higher than the heat losses due to the clas- 
sical mechanism of pair collisions. However, though 
the available theory can predict the conditions 
for development of the drift temperature insta- 
bility, the anomalous heat conductivity coefficient 
can still be calculated only to an order of magnitude 
accuracy. Very roughly, the anomalous heat con- 
ductivity coefficient in Tokamaks can be estimated 
to be 10-20 times as high as the coefficients predicted 
by the classical theory taking into account the role 
played by the confined particles. The experimen- 
tally measured plasma heat losses lie somewhere 
between the classical results and the calculated 
anomalous heat losses. Therefore, the drift instabi- 
lity theory cannot be said to be well substantiated 
by experiments. This inconsistency of the theory 
and experiment reflects the current state of research 
into plasma instabilities. The theoretical classi- 
fication of various instability types in the plasma 
has been in the main completed. The mechanisms 
of instabilities have been analyzed and the condi- 
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tions have been found under which small perturba- 
tions for a given instability decrease or increase. 
However, the linear-approximation analysis of 
such processes is only the first stage of the theory. 

The next, much more difficult stage is to determine 
the coefficients of diffusion and heat conductivity 
in plasma with developing instability. Such a plas- 
ma is typically referred to as the turbulent plasma. 
There has been made only a preliminary analysis of 
diffusion and heat conductivity in the turbulent 
plasma yielding primarily qualitative results which 
are difficult to compare with experimental results. 
The only processes that are comparatively simple 
to analyze are transport phenomena under the 
effect of magnetohydrodynamic instabilities. In 
this case large regions of plasma must move with 
the velocities of the order of the thermal velocity of 
ions. If the magnetohydrodynamic instabilities 
are eliminated, all the processes in the plasma are 
smoothened down. Only drift-type instabilities 
remain, they give rise to plasma turbulence and 
increase the rates of the transport processes. To 
make, at least, rough estimates of the rates of diffu- 
sion and heat transfer daring assumptions have 
to be made, for instance, the hypothesis of the 
"relay" transfer of perturbations in the plasma from 
one magnetic surface with closed field lines to 
another magnetic surface of the same type. Howe- 
ver, at the very best, such an approach yields only 
qualitative results. 

On the other hand, we should note that we still 
cannot make accurate enough measurements of ,*he 
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diffusion and heat conductivity coefficients in 
experiments. To find these parameters we must 
know the distributions of density and temperature 
over the cross section of the plasma and to compare 
the measurements with the predicted values, we 
have to know in detail the geometry of the mag- 
netic field in the plasma. The available experimen- 
tal techniques are largely unsuitable for that. The- 
refore, we have to be content with determining 
such general characteristics as the mean lifetimes 
of particles in the plasma and the mean time of 
energy retention in the plasma. These characteristics 
must be compared with the predicted values proceed- 
ing from fairly arbitrary assumptions about distribu- 
tion of density and temperature in the plasma. So 
the reader can see that the physics of plasma is still 
a rather crude science and it is only the variety of 
sophisticated equations derived by theoreticians 
that to a certain extent can mask its immaturity. 
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